D-A159  729  ATMOSPHERIC  ELECTRIC  I  TV  HAZARDS  THREAT  ENVIRONMENT 

DEFINITIONS)  BOEING  MILITARY  AIRPLANE  CO  SEATTLE  UA 
B  G  MELANDER  ET  AL  AUG  85  AFUAL-TR-85-2052 
UNCLASSIFIED  F33615-82-C-2406  F/G  i/2 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  ST  AN  DARDS  *  1963-A 


OTIC  £iLE  CCPJ 


AFWAL-TR-85-3052 


AD-A159  739 


ATMOSPHERIC  ELECTRICITY  HAZARDS  THREAT 
ENVIRONMENT  DEFINITION 


B.G.  MELANDER 

BOEING  MILITARY  AIRPLANE  COMPANY 
SEATTLE,  WASHINGTON  98124 

AUGUST  1985 

FINAL  REPORT  FOR  PERIOD  AUG  82 -MARCH  1984 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


FLIGHT  DYNAMICS  LABORATORY 
AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OHIO  45433 


DTTC 

ELECTEI 

OCT  04  1985! 


85.  lo  03  135 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related  Government  procurement  operation, 
the  United  States  Government  thereby  incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications ,  or  other  data,  is  net  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture 
use,  or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (ASD/PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS,  it  will 
be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


SAMUEL  E.  HAGINS,  2Lt ,  USAF 

AEHP  Project  Engineer 

Advanced  Development  Programs  Branch 


WALLACE  C.  BUZZARD,  Chief 
Advanced  Development  Programs  Branch 
Vehicle  Equipment  Division 


FOR  THE  COMMANDER 


LAWRENCE  L.  MIDODP  / 

Acting  Chief 

Vehicle  Equipment  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or 
if  the  addressee  is  no  longer  employed  by  your  organization  please  notify  AFWAL/FIEA 
W-PAFB,  OH  45433  to  help  us  maintain  a  current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


UNCLASSIFIED  _ 


security  classification  of  this  page 


1,  REPORT  security  classification 

UNCLASSIFIED 


7».  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3b.  DECLASSIFICATION/ DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


#•.  NAME  OF  PERFORMING  ORGANIZATION 

Eoelng  Military  Airplane  Co. 


6c  ADDRESS  (City.  State  and  ZIP  Cod*) 

Seattle,  Washington  98124 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  Public  release; 
distribution  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AFWAL-TR-8 5-3052 


b.  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 


<ir applicable)  Air  Force  Wright  Aeronautical  Laboratories 

Flight  Dynamics  Laboratory  (AFWAL/FIEA 


7b.  ADDRESS  (City.  State  and  ZIP  Code) 

Wright-Patterson  AFB,  OH  45433 


Ba.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 


18b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

(If  applicable) 


10.  SOURCE  OF  FUNDING  NOS. 


Flight  Dynamics  Laboratory  |  AFWAL/FIEA  j  F33615-82-C-3406 

Be.  ADDRESS  (City.  State  and  ZIP  Code) 

Wright-Patterson  AFB,  OH  45433 


11.  TITLE  ( Include  Security  Classification)  AtHlO Spheric 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO 

63205F 

2507 

02 

01 

ion 

12.  PERSONAL  AUTHOR(S) 

Melander,  B.  G.  »  Cooley,  W.  W. 


13a.  TYPE  OF  REPORT  13b.  TIME  COVERED 

Interim  from  Aug  82 


IB.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Yr.,  Mo..  Day)  15.  PAGE  COUNT 

from  Aug  82  to  Mar  84  August  1985  1^5 


m 

tm 

IK 


COSATI  COOES 


GROUP 


01 


1&  SUBJECT  TERMS  (Continue  on  revert*  if  necessary  and  Identify  by  block  number) 

Aircraft  Lightning  Threat  Atmospheric  Electricity  ^ 


%9.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  document  provides  definition  of  the  atmospheric  electricity  threat  environment  as 
part  of  the  Atmospheric  Electricity  Hazards  Protection  (AEHP)  program  under  Air  Force 
contract  F33615-82-C-3406,  This  issue  contains  an  initial  atmospheric  electricity  threat 
to  be  used  for  aircraft  based  on  measurement  of  (1)  the  lightning  threat  at  the  ground, 
and  (2)  the  static  electrification  threat  at  altitude.  The  experimental  basis  for  this 
threat  is  critically  reviewed.  Recommendations  are  made  for  future  work  in  atmospheric 
electricity  environments  to  refine  and  improve  understanding  of  aircraft  AEH  environments. 


/-  / 


30.  DISTRIBUTION/AVAI LABILITY  OF  ABSTRACT 
UNCLASSIFIED/UNLIMITEO  CB  same  AS  RPT.  □  OTIC  USERS  □ 


23s.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

2Lt  S.  Hagins 


DO  FORM  1473, 83 


31.  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


32b  TELEPHONE  NUMBER  22c.  OFFICE  SYMBOL 

I Include  Area  Code) 

(513)  255-2527  AFWAL/FIEA 


(513)  255-2527 


EDITION  OF  1  JAN  73  IS  OBSOLETE.  UNCLASS TF TBT> _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


FOREWORD 


This  document  was  prepared  under  sponsorship  of  the  Air  Force  Writrhf-AernnauMce] 
Laboratory  of  the  Atmospheric  Electricity  Hazards  Protection  ADP,  Vehicle  Equipment 
Branch,  Wright-Patterson  Air  Force  Base,  Ohio  under  Contract  F33615-82-C-3406. 

Rudy  C.  Beavin,  the  AF  Program  Manager,  provided  industry  and  government  comments 
on  the  threat.  Dr  W.  W  Cooyey  was  the  AEHP  program  Manager. 

B.  G.  Mela nder  of  the  Boeing  Military  Airplane  Company  Mechanical /Electrical 
Technology  group  authored  the  document.  Valuable  technical  discussions  aiding 
In  the  preparation  Involved  the  following  people:  W.  P.  Geren  and  D.  F. 
Strawe  of  Boeing,  R.  L.  Gardner  of  MRC,  M.  A.  Uman  of  University  of  Florida, 
and  E.  P.  Krlder  of  University  of  Arizona.  Special  thanks  to  G.  A.  Dubro  and 
P.  L.  Rustan  for  their  Interest  and  comments.  Thanks  also  to  useful  comments 
made  during  the  Anaheim  Aerospace  Industry/Scientific  Community  Review. 


Accession  For 

NT IS  GRA&I 

dtic  tab 

Utvinnow'-ed 

Justification- 


.'.ion/  _ 

Llity  Codes_ 
•  an-l/or 

nu octal 


Cop, 

l  MtPtcreo 


*•  V-  .S  •  * 


TABLE  OF  CONTENTS 


PAGE 


1.0  INTRODUCTION  1 

1.1  OBJECTIVE  AND  SCOPE  1 

1.2  THREAT  SUMMARY  2 

2.0  REVIEW  OF  AVAILABLE  LIGHTNING  DATA  7 

2.1  INTRODUCTION  7 

2.2  DATA  SOURCES  8 

2.3  DATA  APPLICABILITY  8 

2.3.1  Berger,  Garbagnatl  9 

2.3.2  Uman  and  Krlder  10 

2.3.3  Others  11 

2.4  CONCLUSIONS  12 


3.0  INITIAL  LIGHTNING  CURRENT  THREAT  DEFINITION 

3.1  LIGHTNING  THREAT 

3.2  RATIONALE 

3.2.1  Peak  Rate-of-Rlse 

3.2.2  Peak  Current 

3.2.3  Action  Integral 

3.3  COMPARISON  TO  PRIOR  LIGHTNING  STANDARDS 

3.4  LIGHTNING  THREAT  COMPARISON  TO  MEASURED  DATA 

3.4.1  Uman-Krlder  Measured  Field  Data  Comparison 

3.4.2  Current  Rise  Rate 

3.4.3  Peak  Current 

3.4.4  Action  Integral 

3.4.5  Rise  Time 

3.4.6  Fall  Time  to  Half  Peak 

3.5  LIGHTNING  THREAT  TEST  WAVEFORM 

3.5.1  Recommended  Test  Waveforms 

3.5.2  Comments  on  Uman's  Threat  Waveform 

3.5.3  Waveform  Sensitivity 

3.5.4  Threat  Comparison 

3.5.5  High  Frequency  Content 


24 

24 

25 
25 
25 

25 

26 
26 
27 

27 

28 
29 

29 

30 

30 

31 
31 

31 

32 
32 


4.0  LiarmiNG  MODELING  69 

4.1  TRANSMISSION  LINE  MODELS  70 

4.2  CURRENT  MODELS  DERIVED  FROM  FIELD  DATA  73 

4.3  BASIC  PHYSICS  MODELS  74 

4.3.1  Strawe  74 

4.3.2  Gardner  77 

5.0  STATIC  ELECTRIFICATION  87 

5.1  NOISE  SOURCES  88 

5.2  CHARGING  PARAMETERS  89 

5.3  CORONA  90 

5.4  STREAMERS  91 

5.5  COUPLING  OF  NOISE  TO  AIRCRAFT  SYSTEMS  92 

6.0  PARAMETRIC  THREAT  LEVELS  104 

6.1  LIGHTNING  104 

6.2  STATIC  ELECTRIFICATION  106 

7.0  METEOROLOGICAL  PHENOMENON  114 

7.1  Lightning  Flash  Rate  Density  114 

7.2  MontM”  Variation  114 

?.*  T'inma1.  changes  J16 

7,4  Altitude  Statistics  117 

.7,5  '’'eogra^Mc  Effects  117 

7.6  Other  Parameters  118 

7.7  Strike  Probability  to  Aircraft  119 

8.0  SUMMARY  142 

8.1  CONCLUSIONS  142 

8.2  RECOMMENDATIONS  143 


LIST  OF  FIGURES 


m 


FIGURE  NO. 


3.6 

3.7 

3.8 

3.9 

3.10 

3.11 

3.12 

3.13 

3.14 

3.15 


3.18 

3.19 


TITLE 

Lightning  Flash  Striking  an  Aircraft 

Noise  Generation  Sources 

Stepped-Leader  Initiation  and  Propogatlon 

Return-Stroke  Initiation  and  Propogatlon 

Dart-Leader  and  Subsequent  Return  Stroke 

Multi-Stroke  Threat 

Single  Stroke  Threat  Definition 

Definition  of  Single  Stroke  Threat 

Comparison  of  Electric  Field  Data  with 

Initial  Threat  Definition 

Electric  Field  Frequency  Spectra 

Comparison  of  Initial  Threat  Current 

Rise  Rate  with  Available  Statistical  Data 

Maximum  Current  Rise  Rate 

Current  Amplitude  Distributions 

Current  Amplitude  Distributions 

Current  Amplitude  Distributions 

Energy  Input  Distribution 

Energy  Input  Distribution 

Rise  Time  Distributions 

Rise  Time  Distributions 

New  Data  Suggestive  of  Yery  Fast  Rise  Times 

Fall  Time  to  Half  Peak 

Fall  Time  to  Half  Peak 

Severe  AEH  Threat  Versus  "Test-Like"  Have 

Severe  AEH  Threat  Versus  "Test-Like"  Have 

Uman's  Proposed  Lightning  Current  Test 

Standards  (April,  1982) 

Uman's  Threat  Haveform  -  Severe  First 
Return  Strokes 

Uman's  Threat  Haveform  -  Severe  First 
Return  Strokes 


»**  •**  «  vk** 


LIST  OF  FIGURES  (continued) 


FIGURE  NO*  TITLE  PAGE 

3.24  Oman's  Threat  Waveform  -  Severe  First 

Return  Stokes  60 

3.25  Severe  First  Return  Stroke  Threat 

Comparison  -  Time  Domain  61 

3.26  Severe  First  Return  Stroke  Threat 

Comprlson  -  Frequency  Spectra  62 

3.27  Severe  Subsequent  Stroke  Threat 

Comparison  -  Time  Domain  63 

3.28  Severe  Subsequent  Stroke  Threat 

Comparison  -  Frequency  Spectra  64 

3.29  Single  Discontinuous  Waveform  Fit  to 

Uman's  First  Return  Stroke  Threat  65 

3.30  Frequency  Spectra  Comparison  of  Umans's 

First  Return  Stroke  Threat  and  Discontinuous 
Waveform  Fit  66 

3.31  Double  Exponential  Fit  to  Uman's  First 

Return  Stroke  Threat  67 

3.32  Frequency  Spectra  Comparison  of  Uman's  First 

Return  Stroke  Threat  and  Double  Exponential  68 

4.1  Comparison  of  Bruce-Golde  and  Transmission 

Line  Models  to  Data  83 

4.2  Current  Distribution  for  Model  of  Lin,  et.  al.  84 

4.3  Comparison  of  Lin,  et.  al.  Model  to  Data  85 

4.4  Model  II,  Configuration  1  86 

5.1  Charging  Processes  94 

5.2  Noise  Generation  Sources  95 

5.3  P-Statlc  Charging  96 

5.4  Variation  of  Parlcle  Charge  with  Speed  97 

5.5  Advanced  Aircraft  Effect  Intercepting  Area  98 

5.6  Predicted  Charging  Current  for  Advanced  Aircraft  99 

5.7  Measured  Corona  Noise  Spectra  100 

5.8  Typical  Current  Pulse  Induced  by 

Streamer  Discharge  101 


will 


LIST  OF  FIGURES  (continued) 


FIGURE  NO. 


TITLE 


PAGE 


5.9 

5.1 

6.1 
6.2 

6.3 

6.4 

6.5 

6.6 

7.1 

7.2 


7.3 

7.4 

7.5 

7.6 

7.7 

7.8 

7.9 

7.10 

7.11 

7.12 


Streamer  Noise  Spectra  -  Measured  A  Calculated  102 
Measured  Coronal  and  Streamer  Noise  103 

Multi-Stroke  Threat  108 

Lightning  Current  Waveforms  -  Mean  and  Upper 
and  Lower  1%  Yalues  109 

Attached  Lightning  Frequency  Spectral  Bounds  - 
Current  Amplitude  and  Magnetic  Fields  110 

Nearby  Lightning  Magnetic  Field  - 

Frequency  Spectral  Bounds  111 

Measured  Corona  Noise  Spectra  112 

Streamer  Noise  Spectra  -  Measured  and  Calculated  113 

Thunderstorm  Days  (Isokeraunic  Level)  130 

Cumulative  Frequency  Distributions- 

fluinber  of  Strokes  Per  Flash  131 

Diurnal  Variation  of  Thunderstorm  Activitv  132 

Aircraft  Lightning  Strike  Incidents  as  a 

Function  of  Altitude  1  33 

Altitude  Distribuations  of  VHF  Lightning 

Impulse  Sources  ,  134 

Thunderstorm  Penetrations  and  Lightning 

Statistics  135 

Thunderstorm  Penetrations  and  Lightning 

Statistics  136 

Relationship  of  Lightning  Strikes  to  Relative 

Turbulence  and  Precipitation  Intensities  137 

Temperature  Inside  Storms  and  Probability  of 

Direct  Strike  138 

Turbulence  Level  and  Probability  of  Direct  Strike  139 

Rain  Intensity  and  Probability  of  Direct  Strike  140 

Lightning  Flash  Rate  and  Probability  of 

Direct  Strike  141 


ix 


LIST  OF  TABLES 


TABLE  NO. 

TITLE 

PAGE 

1.1 

AEH  Lightning  Threat  Parameters 

3 

3.1 

Present  Lightning  Threats 

35 

3.2 

Lightning  Parameters  Compared  to  Initial  Threat 

36 

3.3 

Mean  Values  of  Lightning  Parameters  Compared 

to  Initial  Threat 

37 

4.1 

Lightning  Channel  Models 

81 

4.2 

Transmission  Line  MOdel  Comparlosn  to 

Ground  Current  Data 

82 

5.1 

Precipitation  Particle  Parameters 

89 

6.1 

Berger  (Negative  Strokes) 

104 

6.2 

Garbagnatl  (Negative  Strokes) 

105 

6.3 

DtVble  Exponential  Fit  to  Berger’s 

M75)  Statistics 

106 

7.1 

Empirical  Relationships  between  Lightning  Flash 

Density  and  Annual  Thunderdays 

1  ?3 

7.2 

Frequency  of  Occurrence  of  Ground  Flashes 

(Visual/Optical  Methods) 

124 

7.3 

Frequency  of  Occurrence  of  Ground  Flashes 

(Recordings) 

125 

7.4 

Frequency  of  Occurrence  of  Ground  Flashes 

(Power  Transmission- line  Studies) 

126 

7.5 

Ratio  of  Cloud  Flashes  to  Cloud-Ground  Flashes 

127 

7.6 

Lightning  Flash  Rate  and  Probability  or  Direct 

115 

7.7 

Frequencies  of  Thunderstorm  Days  in  Florida 

128 

7.8 

Annual  Variation  in  Occurrences  of  Lightning 

129 

and  Rainfall;  Brisbane 
Approximate  Strike  Probability 


7.9 


130 


1 .0  INTRODUCTION 


The  objectives  of  the  Atmospheric  Electricity  Hazards  Protection  (AEHP) 
program  are  to  develop  design  criteria,  guidelines,  and  qualification  test 
procedures  for  mitigating  any  atmospheric  electricity  vulnerabilities  of 
electronic  equipment  in  future  advanced  aircraft  structures.  The  emphasis  is 
on  Indirect  effects  of  lightning/static  electrification  and  their  Induced 
electrical  transients  In  aircraft  wiring  rather  than  direct  arc/spark  effects. 

Potential  hazards  to  electronics  equipment  within  an  aircraft  encountering 
natural  atmospheric  electricity  depends  strongly  on  the  Interaction  of 
airframe  structure  with  a  lightning  strike  and  Its  associated  electromagnetic 
fields.  The  physical  processes  whereby  electrical  transients  are  Induced  by 
lightning  or  static  electricity  within  an  air  vehicle  wiring  Installation  are 
significantly  different  than  processes  previously  considered  for  structural 
damage  due  to  directly  attached  lightning.  The  most  Important  difference  for 
electronic  effects  Is  the  need  to  consider  rate-of-rlse  parameters  for 
currents  and  fields  because  inductive  and  capacitive  coupling  dominates 
Internal  transients  when  circuitry  Is  Isolated  from  vehicle  structure.  Peak 
values  and  time  duration  of  current  and  fields  are  also  important  In 
determining  voltages  and  currents  within  the  structure.  Pulse  energy  is  also 
Important  In  determining  the  rating  for  transient  protective  devices  that  may 
be  required  to  absorb  energy  from  the  lightning  Induced  transients. 

Physical  understanding  of  the  experimental  data  sources  that  define  lightning 
and  static  electrification  environments  are  Important  In  determining  the 
critical  rate  of  rise  parameters.  These  parameters  are  especially  important 
for  determining  currents  and  electromagnetic  fields  associated  with  aircraft 
Interaction  with  the  AEH  threat  environments. 

1.1  OBJECTIVE  AND  SCOPE 

The  objective  of  this  document  Is  to  provide  definition  of  the  AEHP  program 
threat  environment.  This  threat  document  will  be  revised  periodically  to 
provide  the  "best"  current  definition  of  AEH  threats.  This  Issue  sumarizes 
the  Initial  atmospheric  electricity  threat  to  be  used  for  definition  of  the 
lightning  threat  at  the  ground.  The  experimental  basis  for  this  threat  Is 
critically  reviewed.  Comparisons  with  statistical  data  are  made  as  well  as 


comments  on  the  validity  of  this  data.  Comments  on  on-going  lightning/static 
electrification  data  collection  programs  are  made.  This  data  will  Impact  the 
final  threat  definition.  An  Initial  threat  definition  Is  needed  to  support 
the  Phase  I  environmental  Impact  tests.  A  final  AEH  threat  definition  will 
arise  out  of  additional  studies  of  existing  data  as  well  as  new  data  as  It 
becomes  available. 

1.2  THREAT  SUMMARY 

The  problem  of  determining  the  atmospheric  electricity  (AE)  threat  to  aircraft 
In  flight  has  been  examined  for  many  years.  More  recently  due  to  research 
Into  advanced  composite  material  having  poor  electrical  properties  compared  to 
metals.  Interest  has  Increased  In  better  definition  of  the  AE  threat  to 
aircraft  and  electrical /electronic  equipment.  The  severity  of  the  threat  will 
determine  the  equipment  protection  necessary  for  aircraft  all  weather 
operation. 

Natural  atmospheric  electricity  presents  two  separate  hazards:  lightning  and 
static  discharges  due  to  aircraft  charging.  The  mechanism  of  a  natural  cloud- 
to-ground  lightning  discharge  Is  Illustrated  In  Figure  1.1.  When  sufficient 
charge  accumulates  In  the  lower  part  of  a  cloud  to  cause  an  electric  field 
which  exceeds  the  Ionization  threshold  of  air,  an  electrical  discharge  Is 
Initiated  toward  the  earth.  Because  the  discharge  requires  a  finite  amount  of 
charge  and  time  for  the  channel  resistance  to  lower  to  the  arc  phase,  the 
discharge  proceeds  In  a  sequence  of  steps,  pausing  periodically  to  allow  the 
previous  channel  section  to  become  fully  conducting.  This  mechanism  Is  known 
as  the  stepped  (or  step)  leader  process. 

The  natural  electrical  phenomena  occurring  with  lightning  discharges  vary  In 
number  and  Intensity.  A  statistical  basis  Is  needed  to  define  a  threat 
because  of  this  variability.  A  moderate  threat  level  Is  defined  as  the 
expected  levels  from  a  typical  lightning  flash.  Severe  lightning  Is  defined 
as  a  reasonable  worst-case  level  expected  to  occur  during  the  service  life  of 
an  aircraft.  The  maximum  rate-of-rlse  and  peak  current  values  that  represept 
the  severe  lightning  threat  are  currently  subjects  of  on-going  current 
research  Into  lightning  hazards  to  aircraft  In  flight. 

There  Is  currently  Insufficient  statistical  data  from  efforts  to  measure 
lightning  currents  on  aircraft  In  flight  that  allow  definition  of  a  worst-case 


threat.  Recent  In  flight  data  from  the  NASA  Storm  Hazards  program  has  shown  a 
worst-case  value  of  80  kA/ys  for  low  level  current  from  a  cloud-to-cloud 
strike. 

The  best  available  statistical  data  on  lightning  currents  Is  data  measured  on 
the  ground.  Currently  published  data  on  cloud-to-ground  lightning  currents 
measured  on  the  ground  show  that  180  kA/ys  Is  the  largest  rate  of  rise 
directly  measured.  Recently  obtained  data  on  lightning  currents  Inferred  from 
measured  electric  fields  are  as  high  as  400  kA/ys.  Both  of  these 
measurements  are  subject  to  uncertainties  In  the  measuring  techniques.  These 
limitations  will  be  discussed  In  Chapter  2.  The  best  engineering  judgment  at 
present  Is  that  200  kA/ys  adequately  represents  the  worst-case  expected  for 
an  aircraft  In  flight  and  Is  defined  as  the  severe  level  for  AEH  threat. 
Additional  research  Is  needed  to  resolve  the  difference  between  the  defined 
severe  threat  and  the  distant  field  measurements. 


The  defined  parameters  for  the  AEH  lightning  threats  to  aircraft  are  shown  in 
Table  1-1. 

TABLE  1-1  AEH  LIGHTNING  THREAT  PARAMETERS 


SINGLE 

STROKE 

THREAT 

PEAK 

CURRENT 

MAXIMUM 

RATE  OF 

RISE 

TIME 

TO  HALF 
AMPLITUDE 

ACTION 

INTEGRAL 

Moderate 

(Expected) 

20  kA 

■V 

50  kA/ys 

50ysec 

1.5  x  10*  A2  -sec 

Severe 

200  kA 

200  kA/ys 

50ysec 

1.5  x  106  A2  -sec 

Multiple  Stroke  Events  -  Flash  (Moderate  and  Severe) 


o  Transient 

50-500  ysec 

o  Duration  of 

2  sec 

Duration 

Flash 

o  Interstroke 

10-100  msec 

o  Number  of 

24 

Interval 

Strokes 

A  second  threat  to  aircraft  from  atmospheric  electricity  occurs  due  to 
aircraft  charging  effects.  As  an  aircraft  moves  through  the  air.  It  can 
become  charged.  This  can  result  In  discharges  either  by  streamer! ng,  corona 
or  arcs  and  sparks. 


I 

Figure  1.2  shows  these  discharges  from  various  sections  of  an  aircraft.  Arcs  and 
sparks  can  create  direct  damage  effects  especially  In  fuel  tanks.  Proper 
grounding  and  bonding  can  usually  eliminate  these  problems.  Noise  problems  from 
|  corona  or  streamerlng  phenomena  Is  not  so  easily  eliminated.  These  noise  levels 

must  be  established  as  a  part  of  the  total  atmospheric  threat  levels  to 
aircraft.  The  static  electrification  threat  Is  vehicle  dependent  as  described  In 
Sections  6  and  7. 

The  organization  of  this  document  Is  as  follows.  Lightning  threat  is  discussed 
In  Chapters  2  and  3.  Chapter  2  describes  the  literature  review  and  discusses  the 
choice  of  the  best  data  for  threat  definition.  Chapter  3  outlines  the  ground 
lightning  threats,  both  moderate  and  severe,  and  discusses  their  correlation  with 
available  data.  Chapter  4  discusses  the  applicability  of  lightning  models  as  a 
tool  In  extrapolating  the  lightning  threat  to  aircraft  altitudes.  Chapter  5 
reviews  the  static  electrification  threats  from  streamers  and  corona.  Chapter  6 
lists  parametric  threat  levels  for  both  lightning  and  static  electrification. 
Chapter  7  will  discuss  meterologlcal  phenomenon  and  will  be  completed  for  the 
final  threat  documentation.  Chapter  9  summarizes  the  document  content  and  makes 
recommendations  for  future  work. 
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Figure  I.! 

Lightning  Flash  Striking  an  Aircraft 
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Noise  Generation  Sources 


2.0  REVIEW  OF  AVAILABLE  LIGHTNING  DATA 


2.1  INTRODUCTION 

Lightning  Is  a  transient,  high  current  electric  discharge.  The  most  common 
source  of  lightning  arises  from  the  electric  charge  separation  in  ordinary 
thunderstorm  clouds  (cumulonimbus).  Well  over  half  of  all  lightning 
discharges  occur  within  the  cloud  (Intracloud  discharges).  Cloud- to-ground 
lightning  (sometimes  called  streaked  or  forked  lightning)  has  been  studied 
more  extensively  than  other  forms  of  lightning  because  of  Its  practical 
Interest  (e.g.,  as  the  cause  of  disturbances  In  power  and  communication 
systems,  strikes  to  aircraft  and  the  Ignition  of  forest  fires)  and  because  It 
Is  more  easily  observed  with  optical  Instruments.  Cloud-to-cloud  and 
cloud-to-alr  discharges  are  less  common  than  Intracloud  or  cloud-to-ground 
1 Ightnlng. 

Lightning  strikes  to  aircraft  are  thought  to  Involve  both  cloud- to-ground 
discharges  (at  low  altitudes)  and  Intracloud  events  (at  high  altitudes).  The 
lightning  data  discussed  In  the  present  document  Includes  only  cloud-to-ground 
events.  This  data  will  be  used  to  establish  a  ground  based  threat.  A  later 
version  of  this  document  will  examine  the  airborne  threat  which  will  Include 
cloud  events. 

A  typical  cloud-to-ground  discharge  starts  with  a  preliminary  breakdown  within 
a  cloud  followed  by  a  stepped  leader  Initiating  the  first  stroke  (see  Figure 
2.1).  Leader  steps  are  usually  lps  long,  tens  of  meters  In  length,  with  a 
pause  between  steps  of  50ps.  The  typical  leader  current  Is  the  order  of 
1  kA. 

As  the  leader  nears  an  aircraft.  It  enhances  the  local  fields.  Discharges  off 
the  extremities  of  the  aircraft  are  produced  when  the  field  values  reach  air 
Ionization  levels.  When  the  stepped  leader  connects  to  one  of  the  aircraft 
discharges.  It  becomes  merely  a  part  of  one  step  as  the  stepped  leader 
proceeds  to  ground  or  another  charged  cloud. 

As  the  leader  tip  nears  the  ground,  an  upward  moving  discharge  Is  Initiated  at 
the  ground  (see  Figure  2.2).  The  leader  channel  Is  discharged  as  a  ground 
potential  wave,  the  return  stroke,  propagates  up  the  Ionized  leader  path.  The 
return  stroke  has  a  much  higher  current  than  a  leader,  an  average *of  20  kA 


with  maximum  up  to  200  kA  and  a  rise  time  of  few  microseconds.  The  return 
stroke  may  be  followed  by  a  dart  leader  which  Initiates  a  subsequent  stroke 
(see  Figure  2.3).  Subsequent  strokes  have  faster  rise  times  but  comparable 
rates  of  change  to  return  strokes.  Many  additional  complete  discharges  called 
multiple  strokes  can  take  place.  In  general,  these  subsequent  strokes  have 
smaller  magnitudes.  Establishing  the  size  of  these  threats  Is  the  first  part 
of  establishing  an  AEH  threat. 

2.2  DATA  SOURCES 

The  data  needed  for  establshlng  aircraft  lightning  threat  parameters  are 
measurements  of  direct  lightning  strike  currents  and  EM  fields  on  aircraft  at 
altitude.  Lightning  strike  data  measured  on  aircraft  are  not  yet  well 
understood  because  of  the  limited  number  of  validated  measurements.  Two 
recent  programs  to  gather  more  strike  data  at  altitude  are  the  NASA  Storm 
Hazards  Research  F-106  and  USAF  C-130.  This  data  will  be  used  to  aid  In 
understanding  the  threat  level  at  altitude. 

The  available  data  necessary  to  establish  an  aircraft  lightning  threat 
characterization  Is  limited.  Of  all  the  types  of  lightning  processes  (leader, 
dart  leader,  J  1  K  changes,  preliminary  breakdown,  etc.),  the  most  critical 
processes  to  aircraft  safety  are  thought  to  be  first  return  and  subsequent 
strokes  due  to  high  current  levels,  high  current  rise  rates  and  high  energy 
Inputs  Into  aircraft  systems.  Data  on  these  processes  are  for  the  most  part 
confined  to  ground  measurements  of  electric  and  magnetic  fields,  current 
waveforms,  and  stroke  velocities  (see  Chapter  2  reference  list  for  each  topic). 

Ground  parameters  are  thought  to  be  the  worst  case  situation  for  lightning 
threats.  Current  and  field  amplitudes  and  rate  of  rise  are  thought  to 
decrease  with  altitude.  Hence,  the  most  severe  lightning  threat  to  aircraft 
Is  a  severe  lightning  strike  on  the  ground. 

The  Initial  lightning  threat  outlined  In  this  document  Is  based  on  these 
ground  measurements.  The  final  threat  will  take  Into  account  altitude  effects 
Including  In  flight  data. 

2.3  DATA  APPLICABILITY 

The  present  threat  definition  Is  based  on  ground  based  direct  current 

measurements.  The  data  sets  used  were  taken  from  available  statistical 

studies  on  measured  lightning  current  parameters. 
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Lightning  current  parameters  derived  from  distant  EH  field  measurements  were 
not  used  In  the  Initial  threat  definition.  The  statistics  were  not  yet 
available  when  establishing  the  threat. 

The  errors  associated  with  both  direct  current  measurements  and  currents 
derived  from  field  values  should  be  established.  This  will  allow  a  more 
complete  and  thorough  use  of  the  presently  available  data. 

The  following  sections  describe  the  major  experimental  efforts.  The  author 
date  references  are  for  chronology.  Specific  data  sources  are  Identified  in 
the  chapter  where  the  data  Is  discussed. 

2.3.1  Berger,  Garbagnati 

Berger  (1975)  In  Switzerland  and  Garbagnati  (1982)  In  Italy  have  made  a  series 
of  direct  current  measurements  using  towers  on  mountain  tops.  Berger's  data 
Is  measured  on  a  70  m  tower  on  top  of  Mount  San  Salvadore  which  Is  914  m  above 
sea  level.  Garbagnatl's  measurements  are  from  40  m  towers  on  Italian 

mountains  near  San  Salvadore.  Berger's  latest  data  Is  comprised  of  101 
negative  and  26  positive  first  strokes  and  135  subsequent  strokes.  Garbagnati 
has  103  negative  and  5  positive  first  strokes  and  175  subsequent  strokes. 

Uncertainties  associated  with  both  Berger  and  Garbagnatl's  data  arise  from  the 
fact  that  both  used  towers  on  rocky  mountain  tops.  The  presence  of  the  tower 

may  have  two  effects  on  the  data  taken.  The  presence  of  the  tower  may 

Influence  the  statistics  of  the  lightning  strikes.  For  example,  larger 

amplitude  strokes  are  thought  to  strike  tall  towers.  The  distance  over  which 
a  structure  attracts  a  downward  leader  Is  a  function  of  the  charge  on  the 
leader  which.  In  turn.  Is  related  to  the  amplitude  of  the  current  In  the 
return  stroke  (Golde,  1977).  It  would  thus  follow  that  lightning  strikes  to 
open  ground  should  have  a  greater  proportion  of  lower  currents  while  taller 
structures  might  be  expected  to  be  subjected  to  a  higher  number  of  more 
Intense  currents.  Another  statistical  fluctuation  may  be  In  the  number  of 
positive  strokes  seen.  Evidence  exists  that  positive  flashes  may  Increase 
with  altitude  (Berger,  1975;  Erlkson,  1978).  If  this  trend  holds  true  to 
aircraft  altitudes,  then  any  statistics  taken  near  the  ground  will  not  be  the 
statistics  at  altitude.  The  statistics  of  lightning  strikes  to  aircraft  are 
the  statistics  of  Interest  here.  Any  way  of  measuring  lightning  strike 
statistics  other  than  from  strikes  to  the  specific  aircraft  under  mission 
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conditions  will  not  reflect  the  true  statistics  of  Interest.  Obviously, 
comprotrf ses  on  this  Issue  will  have  to  be  made  and  some  estimates  of  the 
errors  Involved  should  be  established. 

The  presence  of  a  tower  In  direct  current  measurements  also  affects  the 

measurements  themselves.  The  maximum  rise  rate  Is  effected  by  the  tower 

Inductance  and  ground  Impedance  effects.  This  would  tend  to  lower  the 

observed  rise  rates.  The  magnitude  of  this  effect  needs  to  be  established. 
The  early  time  portion  of  the  waveform  may  also  be  effected  due  to  the 
presence  of  the  upward-going  leaders  which  may  lead  to  slower  rise  times  for 
tall  objects  (Clanos  and  Pierce,  1972).  Both  effects  could  Imply  the  current 
rise  rates  measured  on  towers  are  too  low. 

The  leader  effect  should  be  more  pronounced  for  first  return  strokes  than  for 
subsequent  strokes  since  subsequent  strokes  are  not  thought  to  have  long 

upward  propagating  leaders.  Both  Berger  and  Garbagnatl 's  data  show  much 

higher  rates  of  rise  for  subsequent  than  first  return  strokes  while  Uman 

(1973)  and  Weldman  and  Krlder  (1978,  1980)  report  no  difference  In  current 
rates  of  rise  derived  from  fields. 

2.3.2  Uman  and  Krlder 

Uman  and  Krlder  have  spent  many  years  measuring  electric  and  magnetic  fields 
of  various  types  of  lightning  (e.g.  Lin  et  al  1979,  Tiller  et  al  1976,  Uman  et 
al,  1976,  Uman  et  al  1973,  Weldman  et  al  1981).  They  have  found  much  higher 
field  rise  times  (under  a  microsecond)  than  have  been  seen  In  direct  current 
measurements.  The  fast  rise  times  were  seen  In  measurements  over  salt  water. 
Propagation  over  salt  water  does  not  attenuate  the  high  frequencies  as 
severely  as  earth.  The  reported  upper  frequency  limit  of  these  measurements 
Is  20  MHz  (Weldman  et  al  1981)  due  to  wave  action  Influencing  measurements 
above  20  MHz.  This  data  Is  very  useful  In  establishing  validity  of  a  physical 
current  model  from  which  EM  fields  can  be  calculated  and  compared  to  the 

measured  results. 

To  use  this  data  In  establshment  of  a  lightning  current  threat,  the 

uncertainty  associated  with  deriving  currents  from  EM  field  measurements  must 
be  established,  first,  uncertainties  arise  from  the  measurements  due  to 
equipment  limitations,  resonance  effects,  propagation  effects,  etc.  Second 
and  more  serious  are  uncertainties  due  to  the  assumptions  In  the  current 


models  and  the  number  of  parameters  needed  to  fit  the  data.  The  first  type  of 
uncertainty  can  be  readily  quantified.  The  second  needs  to  be  examined.  In 
Uman's  model,  current  Is  assumed  to  propagate  up  the  channel  at  constant 
velocity.  The  channel  Is  taken  to  be  vertical  with  the  Initiation  of  the 

stroke  at  the  ground.  Three  types  of  current  profiles  (uniform  or  leader 

current,  breakdown  current  and  corona  current)  have  been  Incorporated  to  fit 
simultaneous  electric  and  magnetic  fields  at  two  distances.  The  pulse 

velocity  Is  an  Input  to  the  model  and  Is  not  well  known.  Assuming  the 

velocity  to  be  constant  with  height  may  also  add  uncertainties  to  the 
results.  The  corona  current  shape  Is  somewhat  arbitrary  as  stated  by  Lin  et 
al  (1980).  The  non-uniqueness  of  the  current  decomposition  Is  a  major  problem 
and  the  uncertainties  due  to  this  have  not  yet  established.  The  assumption 
that  the  Initiation  point  Is  at  the  ground  and  not  some  distance  above  has 
been  argued  (C.  Baum  private  communication)  to  Introduce  a  factor  of  two  error 
due  to  the  two  return  front  waves  produced  above  ground.  Lastly,  the  effect 
of  the  assumptions  that  the  column  Is  vertical  and  straight  must  be  taken  Into 
account.  Due  to  the  above  assumptions,  error  bars  on  measured  field  data  and 
lightning  current  model  assumptions  must  be  established  before  the  data  can  be 
Incorporated  Into  the  AEHP  lightning  threat  definition. 

2.3.3  Others 

Other  data  sources  are  listed  In  the  Chapter  2  reference  list  (taken  from  a 
review  by  Uman  and  Krlder,  1981).  Only  a  few  names  will  be  specifically 
mentioned  here.  The  recent  review  by  Uman  and  Krlder  (1981)  contains  an 
extensive  review  of  the  lightning  literature. 

Data  sunmarlzed  by  Clanos  and  Pierce  (1972)  was  used  In  establishing  the 
Initial  lightning  threat.  This  data  was  a  compilation  of  work  prior  to  1972 
and  Is  limited  by  the  rise  time  resolution  used  In  the  data  collection.  This 
skews  the  results  toward  longer  rise  times  and  lower  rise  rates. 

Data  collected  and  reviewed  by  Popolansky  (1972)  Include  not  only  results 
obtained  on  tall  chimneys  and  lightning  rods  but  also  the  negative  and 
positive  first  strokes  recorded  by  Berger,  totaling  624  waveforms.  The 
resulting  cumulative  frequency  distribution  curve  produces  a  median  value  of 
28  kA.  Berger  et  al  (1975)  conclude  that  the  median  values  obtained  on  Mount 
San  Salvatore  and  on  tall  chimneys  In  open  country  are  similar.  However,  it 
should  be  pointed  out  that  the  slopes  of  the  curves  of  best  fit  for  these  two 

data  sets  do  not  coincide  completely  (see  Figures  3.9  and  3.10,  next  chapter). 
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Anderson  and  Erlkson  (1978)  In  Sooth  Africa  measured  lightning  currents  on  a 
tower  located  In  open  country*  Unlike  Berger  and  Garbagnatl,  the  measurements 
were  taken  on  relatively  flat  terrain.  Only  a  small  number  of  strokes 
(eleven)  were  analyzed  with  a  maximum  current  rise  rate  of  180  kA/ys  for  a 
subsequent  stroke.  This  Is  a  higher  level  both  In  absolute  value  and  relative 
percentage  than  Berger  and  Garbagnatl 's  measurements.  With  so  few  events, 
however,  definite  comparisons  are  premature. 

Recently  current  waveforms  are  being  measured  In  aircraft  In  flight  by  NASA 
(see  Pitts  1981  and  Pitts  and  Thomas  1981).  To  date  cloud  discharges  dominate 
their  results.  The  1983  program  will  attempt  to  measure  some  cloud-to-ground 
strikes. 

Return  and  subsequent  stroke  velocity  data  have  been  measured  using  luminosity 

data  (e.g.  Orville,  1968;  Boyle  and  Orville,  1976;  Lin  et  al ,  1979;  Hubert  and 

Mouget,  1980;  Jordan  and  Uman,  1980;  Weldman  and  Krlder,  1980).  The  data  has 
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a  wide  variation  of  velocities  ranging  from  2  x  10  m/sec  to  2  x  10 
m/sec.  Typically  the  data  only  allows  calculation  of  an  average  velocity 
found  by  knowing  the  distance  traveled  divided  by  the  elapsed  time.  This  data 
can  be  used  to  compare  lightning  model  results  for  height  variation  of 
velocities.  This  data  can  also  be  used  to  estimate  some  of  the  uncertainty 
present  In  calculating  current  waveforms  from  field  measurements. 

2.4  CONCLUSIONS 

To  correctly  assess  the  ground  lightning  threat,  the  available  data  sets 
(Berger,  Garbagnatl,  Uman,  Eriksson,  Clanos  and  Pierce  and  Popalansky)  must  be 
critically  assessed  as  to  accuracy  and  limitations.  The  earlier  data  (Clanos 
and  Pierce  and  Popalansky)  do  not  reflect  the  recent  fast  rise  times  measured 
and  so  are  biased  to  smaller  rise  rates.  Both  tower  current  data  and  current 
values  derived  from  field  measurements  are  subject  to  uncertainties.  These 
uncertainties  need  to  be  Quantitatively  assessed  before  used  as  a  basis  for 
the  final  current  threat  levels. 

2.5  BIBLIOGRAPHY 
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collected  by  M.  Uman  and  E.  P.  Krlder. 
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Return-Stroke  Initiation  and  Propagation 


He)  Dart  leader  deposits  negative  charge  on  defunct  first-stroke  channel. 
I  He)  Return-stroke  propagates  from  ground  to  cloud.  Scale  of  drawing 
distorted  for  illustrative  purposes.  Adapted  from  Uman  (1971). 
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3.0  INITIAL  LIGHTNING  CURRENT  THREAT  DEFINITION 

The  Initial  AEH  lightning  threat  defined  In  this  chapter  will  be  shown  to 
match  well  with  the  most  recently  published  statistical  lightning  data.  Two 
sets  of  parameters  are  used  to  Identify  a  moderate  or  expected  lightning 
stroke  and  a  severe  or  worst-case  stroke.  The  lightning  threat  Is  not 
designed  to  match  the  physical  parameters  of  a  particular  lightning  stroke  but 
rather  to  be  representative  of  the  range  of  values  for  the  many  types  of 
lightning  discharges. 

3.1  LIGHTNING  THREAT 

The  parameters  of  most  Importance  In  this  lightning  threat  were  chosen  because 
of  their  Impact  on  aircraft  electrical/electronic  systems.  These  parameters 
Include  maximum  current  rise  rate,  peak  current  and  energy  Input  (action 
Integral).  The  selected  values  were  chosen  from  a  review  of  existing  data. 
They  are  shown  to  be  consistent  with  statistical  variations  of  other  available 
data  In  the  following  sections.  Other  characteristics  of  the  waveform  (rise 
time  and  fall  time)  are  determined  uniquely  by  choosing  the  above  three 
parameters  since  the  threat  model  has  only  three  Independent  parameters.  Even 
so,  both  rise  and  fall  times  are  well  within  the  statistics  of  measured  data. 

The  lightning  threat  must  characterize  both  single  and  multiple  stroke 
phenomena.  Multiple  strokes  will  be  characterized  by  several  single  strokes 
with  the  addition  of  Induced  transient  duration.  Inter-stroke  time  Interval, 
total  event  time  and  total  number  of  strokes.  These  parameters  are  listed  In 
Figure  3.1. 

The  Initial  single  stroke  lightning  threat  model  is  a  double  exponential 
waveform  representing  the  lightning  current;  the  waveform  and  spectrum  is 
shown  In  Figure  3.2.  The  double  exponential  form  of  the  lightning  threat 
model  will  be  shown  to  adequately  predict  the  expected  electric  field  spectra 
from  a  combination  of  measured  lightning  discharges.  This  threat  Is  to  be 
Interpreted  as  the  current  flowing  In  the  unperturbed  lightning  arc  channel 
(l.e.,  no  aircraft  Interaction).  In  applying  the  threat  to  an  aircraft,  an 
electromagnetic  coupling  model  Is  needed  which  Includes  the  aircraft  and 
channel  geometry  and  Includes  the  threat  as  an  Incident  current  waveform 
propagating  along  the  channel. 


r: 
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The  double  exponential  waveform,  which  Is  a  convenient  model  for  engineering 
calculations,  has  a  historical  precedent  In  both  lightning  and  NEMP  (Reference 
3-1).  The  waveform  parameters  are  selected  to  Include  the  most  Important 
features  of  the  lightning  current  rather  than  to  faithfully  represent  any 
single  lightning  waveform.  The  three  lightning  parameters  Incorporated  In  the 
threat  waveform  are  peak  rate-of-rlse,  peak  current,  and  action  Integral.  The 
moderate  and  severe  threats  were  selected  to  be  expected  and  worst-case  levels 
of  the  parameters  based  upon  review  of  the  best  available  measured  data. 
Figure  3.3  shows  these  threat  parameter  values  and  the  defining  equations  for 
the  double  exponential  current  waveform. 

3.2  RATIONALE 

No  single  waveform  can  re§rtsent  all  types  of  lightning  discharges  (e.g., 
cloud-to-ground.  Intracloud,  positive  strokes,  negative  first  strokes,  and 
negative  subsequent  strokes).  It  Is  necessary  to  select  parameters  from 
particular  stroke  types  which  provide  reasonably  conservative  threat  levels 
for  all  strokes.  Cloud-to-ground  strokes  were  chosen  because  they  are 
generally  more  severe,  although  a  less  frequent  threat  to  aircraft,  than 
Intracloud  discharges.  On  this  basis,  the  threat  waveform  parameter  values 
were  chosen  as  described  below. 

3.2.1  Peak  Rate-of-Rlse 

Rise  rate  data  displayed  by  Clanos  and  Pierce  has  recently  been  Interpreted  to 
be  too  low.  More  recent  statistical  data  [References  3-4,  3-6,  3-8,  and  3-13] 
has  shown  higher  current  rates  of  rise.  The  values  for  the  Initial  threat 
determination  were  taken  from  Berger's  [Reference  3-6]  tower  measurements. 
The  moderate  value  of  50  kA/ys  lie  at  his  upper  35X  mark  while  the  severe 
threat  of  200  kA/ys  was  chosen  at  his  1%  level. 

3.2.2  Peak  Current 

The  peak  current  was  chosen  from  the  statistical  study  done  by  Clanos  and 
Pierce  [Reference  3-7].  The  moderate  threat  level  of  20  kA  was  chosen  at 
their  50%  level  for  first  return  strokes.  The  severe  threat  of  200  kA  was 
chosen  at  their  upper  IX  level. 

3.2.3  Action  Integral 

The  action  Integral  was  also  chosen  from  Clanos  and  Pierce  data.  Peak  current 

values  and  mean  rise  and  fall  times  were  used  to  determine  a  moderate  energy 

Input  level  of  1.5  x  10*  A2-s  and  a  severe  level  of  1.5  x  10®A2-s. 
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3.3  COMPARISON  TO  PRIOR  LIGHTNING  STANDARDS 

In  this  section  the  lightning  threat  defined  In  Figures  3.1  through  3.3  is 
compared  against  other  Industry  lightning  standards.  Table  3.1  summarizes 
comparisons  between  previously  used  standards  and  the  present  threat 
definition.  Only  the  severe  threat  Is  compared  since  the  other  standards 
represent  severe  strikes. 

The  maximum  current  rise  rate  of  Industrial  lightning  threat  standards  are 
given  In  Table  3.1.  The  values  In  parenthesis  are  obtained  from  the  peak 
current  value  divided  by  the  rise  time  for  the  corresponding  standard.  This 
gives  an  average  rise  time  or  a  maximum  rise  time  If  the  Initial  current 
waveform  Is  a  straight  line.  The  straight  line  waveform  for  the  initial 
current  rise  Is  used  In  the  SAE-4L  standards.  The  rise  time  for  the  AEH 
threat  Is  the  peak  rate  of  rise,  not  an  average.  For  a  double  exponential 
waveform  these  quantities  are  substantially  different.  The  peak  rate  of  rise, 
not  the  average.  Is  the  Important  quantity  when  specifying  a  lightning 
threat.  This  difference  between  peak  and  average  rate  of  rise  leads  to  the 
AEH  higher  standard  for  rise  rate. 

The  peak  and  average  rate  of  rise  definition  also  leads  to  an  AEH  rise  time 
different  from  other  Industry  standards.  These  rise  time  differences  are 
shown  In  Table  3.1.  The  longer  AEH  rise  time  is  due  to  the  double  exponential 
waveform  used.  Other  standards  use  linear  ramp  functions.  However,  the  rise 
rate,  not  the  rise  time.  Is  the  Important  parameter  when  considering  possible 
damage  to  electronic  equipment. 

The  peak  current  and  fall  time  are  listed  In  Table  3.1.  The  AEH  initial 
threat  values  agree  with  those  previously  used. 

3.4  LIGHTNING  THREAT  COMPARISON  TO  MEASURED  DATA 

The  data  presented,  for  comparison  between  the  AEHP  defined  threat  and 
lightning  measurements.  Is  ^he  most  recent  found  in  the  literature.  Each 
model  parameter  Is  examined  separately  below.  Two  summary  tables  of 
experimental  statistics  are  given  In  Tables  3.2  and  3.3.  The  first  table  Is 
taken  from  Berger  (Reference  3-2).  The  second  Is  from  Garbagnatl  and  Loplparo 
(Reference  3-4).  Tables  3.2  and  3.3  summarize  the  data  from  References  3-2 
and  3-4  from  which  most  of  the  data  comparisons  In  this  section  were  made. 


All  the  AEHP  threat  parameters  are  within  reasonable  statistics  of  this 
database.  Detailed  comparisons  with  recent  Individual  experlmentors  are  given 
In  the  following  subsections. 

3.4.1  Uman-Krlder  Measured  Field  Data  Comparison 

A  comparison  Is  made  In  this  section  to  measured  electric  fields  presented  by 
Uman  and  Krlder  (Reference  3-7).  The  comparison  of  a  lightning  current 
profile  to  electric  radiation  field  data  Is  not  straightforward  If  lightning 
geometry,  propagation  effects,  height  dependence  of  lightning  current,  etc. 
are  taken  Into  account.  Since  the  comparison  presented  here  Is  meant  to  be 
preliminary,  a  simpler  procedure  will  be  followed.  The  magnetic  radiation 
field  Is  calculated  from  an  Infinite  current  carrying  wire.  This 

approximation  Is  valid  close  to  the  wire  (less  than  a  wavelength  away).  This 
minimum  distance  Is  300  km  fori  khz  and  goes  to  3  Km  at  loo  khz.  The  electric 
radiation  field  Is  obtained  from  the  magnetic  radiation  field  assuming  free 
space  radiation.  The  results  for  both  moderate  and  severe  current  profiles 
are  plotted  In  Figure  3.4  against  Uman  and  Krlder's  (Reference  3-7)  electric 
field  spectra  for  first  strokes.  Also  plotted  are  NEMP  results  for 

comparison.  The  threat  models  tend  to  be  well  above  the  data  for  low 
frequency  and  straddling  the  data  at  higher  frequency.  The  shape  of  the 
frequency  spectra  of  the  predicted  fields  versus  the  measured  fields  Is  fairly 
good.  Of  course  this  assumes  no  dispersion  effects  occur  In  propagation. 
This  Is  not  the  case  especially  at  higher  frequencies  or  where  the  earth  is 
very  lossy.  Newer  measurements  by  Krlder  are  over  salt  water  having  low  loss 
propagation. 

Figure  3.5  taken  from  Reference  3-7  shows  frequency  spectra  from  different 
lightning  discharges.  The  solid  line  represents  the  first  return  stroke 
data.  All  the  discharges  shown  are  on  top  of  or  close  to  the  first  return 

stroke  data.  This  justifies  our  use  of  AEH  threat  comparison  to  first  stroke 

data  In  Figure  3.4. 

3.4.2  Current  Rise  Rate 

The  single  stroke  threat  model  value  for  maximum  current  rise  rate  Is  given  In 
Figure  3.3  as  50  kA/ys  for  a  moderate  stroke  and  200  kA/ys  for  a  severe 
one.  The  values  correspond  to  the  upper  35%  and  upper  1%  as  shown  In  data 
from  Berger  (Reference  3-2)  In  Figure  3.6.  Note  that  no  measurements  were 
made  with  rate  of  rise  greater  than  100  kA/ys.  Figure  3.7  shows  older  data 
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accumulated  by  Clanos  and  Pierce  (Reference  3-5).  The  moderate  and  severe 
threat  correspond  to  the  upper  12%  and  upper  0.1%  values.  The  Clanos  and 
Pierce  data  Is  actually  average  rates  of  rise  rather  than  maximum  values. 
Thus  the  threat  values  show  as  higher  percentiles  as  compared  to  Berger's 
data.  Data  from  Garbagnatl  (Reference  3-4)  of  two  different  data  sets  (first 
and  subsequent  strokes)  Is  shown  In  Figure  3.8.  The  moderate  threat  Is  at  the 
10%  mark  In  (a)  and  35%  In  (b).  The  severe  threat  is  at  3%  In  (a)  and  1.5%  In 
(b).  These  statistical  values  correlate  well  with  the  choice  of  rise  rate 
picked  In  the  current  threat  model. 

There  Is  considerable  uncertainty  In  the  maximum  rate  of  rise  values.  It 
should  be  noted  that  the  largest  current  rise  rates  measured  by  Berger  or 
Garbagnatl  and  Loplparo  were  102  kA/ysec.  This  data  may  be  low  because  of 
Instrumentation  bandwidth  limitations  and  the  tower  heights.  No  data  has  been 
directly  measured  showing  current  rise  rate  greater  than  the  180  kA/ys 
reported  by  Anderson  and  Erickson.  Anderson  and  Erlkson  (1978)  In  South 
Africa  measured  lightning  currents  on  a  tower  located  In  open  country.  Unlike 
Berger  and  Garbagnatl,  the  measurements  were  taken  on  relatively  flat 
terrain.  Only  a  small  number  of  strokes  (eleven)  were  analyzed  with  a  maximum 
current  rise  rate  of  180  kA/ys  for  a  subsequent  stroke.  This  Is  a  higher 
level  both  In  absolute  value  and  relative  percentage  than  Berger  and 
Garbagnatl's  measurements.  With  so  few  events,  however,  definite  comparisons 
are  premature. 

Higher  rate  of  rise  values  (up  to  400  kA/ys)  have  been  reported  (Reference 
3-8)  but  are  Inferred  from  distant  electric  fields. 

Consequently,  the  AEHP  threat  value  has  been  adjusted  upward  to  200  kA/us. 
This  Is  based  on  engineering  judgement  of  where  the  real  threat  is  expected  to 
be.  Further  testing  Is  needed  to  resolve  the  threat  uncertainty  but  It  is 
expected  to  be  within  a  factor  of  two  of  the  present  AEHP  threat. 

3.4.3  Peak  Current 

The  peak  current  values  were  taken  from  Poplansky's  data  (Reference  3)  shown 
In  Figure  3.9.  The  moderate  (20  kA)  and  severe  (200  kA)  threat  values  were 
chosen  at  60%  and  1%  respectively.  Data  In  Figure  3.10  from  Berger  (Reference 
3-2)  show  these  values  to  lie  at  80%  and  0.1%,  respectively.  Figure  3.11  Is 

data  from  Garbagnatl  (Reference  3-4).  The  moderate  threat  Is  at  80%  In  9(a) 
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and  60%  In  9(b).  The  severe  threat  Is  less  than  0.1%  In  both  figures.  The 
moderate  threat  Is  close  to  a  median  value  In  Figures  3.9  and  3.11  and  seems 
to  be  a  good  choice.  All  the  data  but  Poplansky's  have  smaller  current  values 
near  the  severe  threat  limit.  Our  use  of  200  kA  seems  justified,  at  least 
Initially,  on  the  basis  of  prior  Industry  standards  and  because  the  200  kA  is 
above  the  1%  levels  of  Figures  3.9  and  3.11  and  so  Includes  at  least  the  1% 
level  In  the  other  data  sets  available. 

3.4.4  Action  Integral 

The  last  parameter  chosen  in  the  threat  definition  Is  the  energy  input  or 
action  Integral.  Figure  3.12  shows  the  Berger  data.  The  percentiles  for  the 
moderate  threat  energy  case  range  from  25%  to  85%  for  negative  subsequent  and 
first  return  strokes.  The  severe  threat  Is  at  1%  of  the  first  negative  stroke 
and  <  0.1%  for  the  subsequent  stroke.  The  severe  threat  level  agrees  with 
the  first  return  stroke  and  Is  overly  severe  for  subsequent. 

Figure  3.13  presents  data  from  Reference  3-4.  3.13(a)  has  moderate  and  severe 
threat  energy  levels  at  82%  and  1%  while  3.13(b)  has  48%  and  <  0.1%  levels. 
The  moderate  threat  agrees  with  the  subsequent  stroke  data  and  is  not  severe 
enough  according  to  first  return  stroke  data.  The  severe  threat  level  is 
overly  severe  for  the  subsequent  and  adequate  for  first  return  strokes. 

3.4.5  Rise  Time 

The  rise  time  of  the  threat  Is  a  consequence  of  fixing  the  maximum  current 
rise-rate,  the  peak  current  and  the  action  Integral.  Rise  time  Is  not  fixed 
Independently  of  these  parameters.  Consequently,  comparing  rise  time  to  data 
gives  an  Indication  of  the  general  adequacy  of  the  form  of  the  double 
exponential  used  to  model  the  lightning  current.  Figure  3.14  from  Reference 
3-2  shows  the  moderate  threat  value  of  2ys  falling  at  25%  and  95%  for 
subsequent  and  first  strokes  respectively.  The  severe  threat  of  4ys  lies  at 
the  6%  and  70%  mark.  Figures  3.15a  and  b  show  results  from  Reference  3-4. 
Figure  3.15a  shows  moderate  and  severe  percentile  of  70%  and  43%.  Figure 
3.15b  shows  values  of  25%  and  3%.  The  data  Is  widely  scattered  over  a  large 
range  of  rise  times.  The  values  chosen  fall  within  this  middle  range  of  the 
data. 

Figure  3.16  from  Reference  3-6  shows  results  that  have  very  short  rise  times 
(<lys).  These  results  were  taken  over  salt  water.  The  relationship  to 
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|  data  taken  over  land  Is  not  completely  understood  but  the  rise  time  should  be 

slower  over  land  than  over  water.  This  data  Indicates  that  much  shorter  rise 
times  may  have  to  be  Included  In  a  future  threat.  The  Impact  of  short  rise 
I  times  on  the  lightning  threat  Is  through  the  related  parameter  rise  rate. 

Generally  the  shorter  the  rise  time,  the  higher  the  rise  rate,  although  the 
relationship  Is  waveform  dependent.  Current  rise  rate  values  directly  affect 
the  amount  of  protection  an  aircraft  must  have  to  guard  against  lightning 
strike  cause  equipment  malfunctions. 

3.4.6  Fall  Time  to  Half  Peak 

The  fall  time,  like  the  rise  time,  Is  determined  by  the  preceedlng  parameters 
so  the  fall  time  was  not  fit  to  the  measured  data  In  defining  the  threats. 
Figure  3.17  from  Reference  3-2  shows  the  moderate  and  severe  threat  time  of 
50ys  to  lie  at  30%  and  80%  for  negative  subsequent  and  negative  first 
strokes  respectively.  Figure  3.18  has  data  taken  from  Reference  3-4.  The 
AEHP  threat  fall  time  ranges  from  the  35%  to  75%  level  In  15(a)  and  Is  30%  In 
15(b).  The  threat  value  Is  well  within  the  median  range  of  the  above  data. 

3.5  LIGHTNING  THREAT  TEST  WAVEFORM 

Uman  (Reference  3-8)  proposed  lightning  threat  waveforms  for  both  first  return 
and  subsequent  strokes.  A  separate  waveform  for  testing  must  be  used  since 
the  double  exponential  form  cannot  easily  be  produced  In  a  laboratory.  The 
double  exponential  peak  rate  of  rise  may  be  limited  by  stray  Inductance  and 
capacitance  In  a  test  configuration  since  the  peak  occurs  at  zero  time.  The 
purpose  of  this  section  Is  to  examine  the  sensitivity  of  the  time  domain 
current  waveform  on  the  frequency  spectral  components.  The  second  aim  is  to 
compare  the  Initial  AEHP  lightning  threat  waveforms  with  both  Uman's  proposed 
threats  and  a  recommended  test  waveform. 

The  choice  of  a  representative  waveform  of  a  lightning  threat  useable  for  test 
and  analysis  must  consider  both  time  and  frequency  domain  profiles.  The  use 
of  a  piecewise  continuous  waveform  can  Introduce  high  frequency  nulls  which 
are  undeslreable  for  most  analyses  because  of  the  lack  of  frequency  content 
and  hence  response.  The  use  of  a  smooth  contl nous  function  for  analysis  Is 
desired  to  generate  a  smooth  continuous  frequency  spectrum  with  no  such 
nulls.  Real  test  waveforms  will  be  smoothed  by  the  generator  turn-on  and 
parasitic  R/L/C  elements.  Double  exponentials  have  been  shown  to  bound  In 


30 


both  time  and  frequency  domains  more  complex  time  domain  waveforms.  The 
purpose  of  a  lightning  test  waveform  Is  to  bound  both  In  time  and  frequency 
the  Important  parameters  of  lightning. 

3.5.1  Recommended  Test  Waveform 

The  recommended  test  waveform  shape  Is  a  sine  wave  with  an  exponential  tall. 
A  comparison  of  the  recommended  test  waveform  with  the  double  exponential 
threat  waveform  Is  seen  In  Figures  3.19  and  3.20.  The  time  domain  waveforms 
are  shown  In  the  first  figure.  The  comparisons  are  done  for  a  severe  threat. 
The  frequency  spectra  Is  shown  In  Figure  3.20.  The  double  exponential  and 
test  waveforms  match  well  In  the  frequency  domain.  Thus,  the  proposed  test 
waveform  should  represent  the  double  exponential  threat  well. 

3.5.2  Comments  on  Uman's  Threat  Waveform 

Figure  3.21  shows  Uman's  proposed  lightning  current  waveforms  for  both  first 
return  and  subsequent  stroke  (Reference  3-8).  The  waveforms  both  go  to  zero 
at  t  *  0  and  t  »  300  ys.  The  maximum  rise  rates  for  the  two  waveforms  are 
750  kA/ys  for  the  severe  first  return  stroke  and  600  kA/ys  for  the  severe 
subsequent  stroke.  Uman's  severe^threat  rise  rate  values  were  chosen  to  be  a 
factor  of  five  larger  than  the  moderate  rise  rate  values  derived  from  field 
data. 

3.5.3  Waveform  Sensitivity 

The  sensitivity  of  spectral  content  was  examined  for  various  choices  of  time 
domain  waveforms  similar  to  Uman's  threat.  Uman's  severe  first  return  stroke 
profile  Is  shown  In  Figure  3.22  (solid  line).  The  first  change  to  this 
waveform  was  to  delete  the  ramps  running  from  t  ■  0  to  t  *  100  ys  and  from 
t  *  200  to  t  *  300  ys  and  then  to  normalize  the  current  amplitude  to  zero  at 
t  *  100  and  t  *  200  ys.  The  maximum  rise  rate  was  left  unchanged.  The 
second  alteration  was  to  fill  In  the  notch  present  after  the  maximum  amplitude 
Is  reached.  This  change  Is  shown  by  a  dotted  line  near  the  peak  amplitude  In 
Figure  3.22.  Again  the  maximum  rise  rate  was  left  unchanged.  The  spectra  of 
the  waveforms  are  shown  In  Figures  3.23  and  3.24.  Figure  3.23  shows  a 
comparison  of  the  original  waveform  to  one  without  the  ramps.  The  only  change 
Is  at  the  low  frequency  end  of  the  spectrum  as  expected.  No  difference  above 
10  KHz  exists.  Figure  3.24  shows  the  waveforms  with  and  without  the  notch 
present.  Again  there  Is  very  little  difference  over  the  entire  spectrum. 
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3.5.4  Threat  Comparison 

Comparison  of  the  Initial  AEHP  threat  with  Uman's  waveforms  and  a  sample  of 
Berger's  (1975)  data  selected  by  Uman  Is  shown  In  Figure  3.25-3.28.  Figures 

3.25  and  3.27  show  the  time  domain  waveforms  while  Figures  3.26  and  3.26  show 
the  associated  spectra.  Both  the  Initial  AEHP  threat  and  Berger's  data 
envelope  Uman's  threat  waveforms  except  at  early  time.  Uman  picked  750 
kA/ys  for  his  first  return  stroke  threat  and  600  for  his  subsequent  stroke 
threat.  Berger's  example  has  values  of  70  kA/y  for  first  return  stroke 
(upper  20%  of  his  data)  and  360  kA/ys  for  subsequent  strokes  (extremely 
severe,  less  than  upper  0.1%  of  his  data).  The  AEHP  threat  values  are  200 
kA/ys  for  severe  threats  and  50  kA/ys  for  moderate  threats  based  on  both 
first  return  and  subsequent  stroke  data.  Upper  bounds  and  reliable  statistics 
for  high  rates  of  rise  are  not  yet  well  established.  Therefore,  Uman's  high 
rise  rates  may  well  be  Indicative  of  a  more  severe  lightning  stroke  than  the 
upper  1%  level  taken  for  the  basis  of  the  AEHP  threat. 

These  high  rise  rates  are  also  reflected  In  the  frequency  spectra  In  Figures 

3.26  and  3.28.  The  Initial  AEHP  threat  bounds  both  Berger's  data  and  Uman's 
waveform  at  lower  frequencies  for  both  first  return  and  subsequent  strokes. 
At  higher  frequencies,  Uman's  threat  Is  higher  than  either  Berger's  data  or 
the  AEHP  threat  level.  Again  this  Is  due  to  the  rise  rates  chosen. 

Uman's  large  rise  rate  values  are  chosen  from  current  parameters  derived  from 
EM  field  measurements.  These  values  are  subject  to  uncertainty  due  to  the 
deconvolution  from  fields  to  currents.  This  process  Is  nonunique  and  many 
assumptions  must  be  made  as  to  the  current  waveforms.  The  lower  AEHP  threat 
levels  are  taken  from  tower  measurements.  These  values  are  also  subject  to 
uncertainty  due  to  tower-lightning  Interactions  as  discussed  In  the  last 
chapter.  Quantitative  evaluations  of  both  Uman's  field  data  and  the  tower 
measurements  must  be  done  to  establish  a  better  statistical  base  for  rise  rate 
val ues. 


3.5.5  High  Frequency  Content 

The  last  point  to  be  addressed  here  regards  the  high  frequency  shape  of  Uman's 
waveforms.  Above  10  Wz  a  sin  x/x  type  variation  appears  In  the  spectra  with 
deep  nulls  present.  This  type  of  waveform  Is  not  useful  for  either  testing  or 
analysis  due  to  these  nulls.  Rather  what  Is  needed  are  bounding  curves  In  the 
time  and  frequency  regime. 
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The  sin  x/x  variation  Is  caused  by  the  multiple  discontinuities  In  slope  In 
Uman's  time  domain  waveforms.  For  more  than  one  discontinuity  In  slope*  the 
sin  x/x  variation  appears  with  the  first  frequency  null  determined  by  the  time 
spacing  between  the  discontinuities.  An  example  of  this  Is  the  waveform  given 
In  Figure  3.29  shown  again  Uman's  severe  first  return  stroke  threat.  The  two 
frequency  spectra  are  plotted  against  each  other  In  Figure  3.30.  The 
discontinuous  waveform  plotted  has  nulls  Identical  to  Uman's  waveforms  because 
the  rise  times  of  the  steepest  slope  are  Identical  (both  waveforms  steepest 
slope  rise  Is  0.1  ps  -  this  gives  the  first  null  at  10  MHz). 

Rather  than  smoothing  the  discontinuities  In  Uman's  waveforms  which  Is 
tedious,  a  better  approach  Is  to  bound  the  waveform  and  spectrum.  An  example 
of  such  a  bound  by  a  double  exponential  waveform  (but  not  the  Initial  AEHP 
threat  waveform)  Is  given  In  Figures  3.31  and  3.32.  The  parameters  chosen  for 
this  double  exponential  were  chosen  to  bound  Uman's  waveform  at  high 
frequencies.  The  general  expression  with  the  values  chosen  are  shown  below. 

I(t)  ■  I0  [e"01*  -  e"0t]  (1) 

IQ  «  180  kA;  a  *  1.4  E  ♦  4  sec"1;  and  g  «  9.0  E  +  6  sec"1.  (2) 

The  peak  amplitude  of  this  waveform  Is  178  kA.  The  maximum  rate  of  rise  Is 
very  high  at  1600  kA/ps.  The  spectral  fit  to  Uman's  waveform  Is  excellent 
at  high  and  low  frequencies.  Between  10  KHz  and  1  MHz  there  Is  up  to  a  factor 
of  three  difference  In  amplitude.  This  double  exponential  serves  as  an 
excellent  example  of  a  bounding  wave  in  both  time  and  frequency  domains. 


CHAPTER  3  REFERENCES 


Golde,  R.H.,  ed;  "Lightning,  Vol  Is  Physics  of  Lightning",  Academic 
Press,  1977. 

Berger,  K;  Anderson,  R.B.;  Kronlnger,  H.,  "Parameters  of  Lightning 
Flashes",  Electra,  80,  23-37,  1975. 

Popolansky,  F.,  "Frequency  Distribution  of  Amplitudes  of  Lightning 
Currents",  Electra,  22,  139-147,  1972. 

Garbagnatl,  E.;  Loplpard,  G.B.;  "Lightning  Parameters  -  Results  of  10 
Years  of  Systematic  Investigation  In  Italy",  Proceedings  of 
International  Conference  on  Lightning  and  Static  Electricity;  Oxford, 
England,  March,  1982. 

Clanos,  N.;  Pierce,  E.T.,  "A  Ground  Lightning  Environment  for 
Engineering  Usage",  Stanford  Research  Instltude  Technical  Report;  No. 
1,  Project  1834,  August  1972. 

Clifford,  D.W.;  Krlder,  E.P.;  Uman,  M.A.;  "A  Case  for  Submicrosecond 
Rise-Time  Lightning  Current  Pulses  for  Use  In  Aircraft  Induced-Coupling 
Studies",  IEEE  Symposium  on  Electromagnetic  Compatablllty,  San  Diego, 
CA,  October  1979. 

Uman,  M.A.;  Krlder,  E.P.;  "A  Review  of  Natural  Lightning:  Experimental 
Data  and  Modeling",  IEEE  EMC  Transactions;  p.  79-111,  May  1982. 

Weldman,  C.  D.;  Krlder;  "Submicrosecond  Rlsetlmes  In  Lightning 
Return-Stroke  Fields",  GRL,  F,  955-958,  November,  1980. 


Table  3.1 

Present  Lightning  Threats 


Table  3.2 

lightning  Parameters  Compared  to  Initial  Threat 
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Table  3.3 

Mean  Values  of  Lightning  Parameters  Compared  to  Initial  Threat 
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(Non-Damage  Effects) 


MODERATE  THREAT  IEVERE  THREAT 


Figure  3.2 

Single  Stroke  Threat  Definition 
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Figure  3.3 

Definition  of  Single  Stroke  Threat 
(Damage  Effects)  _ 
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Figure  3.4 

Comparison  of  Electric  Field  Data  with 
Initial  Threat  Definitions 


Figure  3.8 

Maximum  Current  Rise  Rate 


Figure  3.9 

Current  Amplitude  Distributions 


Figure  3.10 

Current  Amplitude  Distributions 
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Figure  3.12 

Energy  Input  Distribution 
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Figure  3.13 

Energy  Input  Distributions 
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Figure  3.14 

Rise  Time  Distributions 


Figure  3.15 

Rise  Time  Distributions 


FIRST  RETURN 


Figure  3.16 

New  Data  Suggestive  of  Very  Fast  Rise  Times  (<1  ns) 
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Figure  3.17 

Fall  Time  to  Half  Peak 


Figure  3.18 

Fall  Time  to  Half  Peak 
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ersus  Test-Like  Wave 


Figure  3.21 

’s  Proposed  Lightning  Current  Test  Standards 

April  1 982 


Figure  3.22 

’s  Threat  Waveform 
First  Return  Strokes 


With  and  Without  Ramp 


Figure  3.23 

Uman’s  Threat  Waveform 
Severe  First  Return  Strokes 


With  and  Without  Notch 


Threat  Waveform 
rst  Return  Strokes 


Figure  3.25 

Severe  First  Return  Stroke  Threat  Comparison 

—  Time  Domain  — 


Figure  3.29 

Simple  Discontinuous  Waveform  Fit  to  Uman’s 

First  Return  Stroke  Threat 


Figure  3.31 

Double  Exponential  Fit  to  Uman’s 
First  Return  Stroke  Threat 


4.0  LIGHTNING  MODELING 

A  model  can  be  defined  as  a  physical  or  mathematical  construct  which 
approximates  to  various  degrees  certain  observed  aspects  of  natural  or  man 
made  phenomena.  In  this  section  we  will  be  concerned  with  a  evaluation  of 
various  mathematical  models  of  lightning  processes.  The  end  product  required 
of  the  lightning  modeling  will  be  the  prediction  of  both  lightning  arc  current 
and  EM  field  values  produced  by  these  processes  at  aircraft  altitudes.  This 
model  Is  needed  because  of  the  spare  data  on  the  statistical  variations  of 
lightning  currents  at  altitude.  Furthermore  this  data  will  require 
considerable  time  to  obtain.  An  analytical  model  will  supplement  the  aircraft 
Inflight  measurements  oy  using  ground  statistics  once  the  extrapolation  to 
altitude  Is  established  using  a  few  measurements. 

There  are  basically  three  levels  of  detail  In  current  mathematical  models  for 
lightning  phenomena: 

o  Physical  analog  model  Including  arc,  clouds,  and  electromagnetic 
propagation 

o  Lumped  constant  electrical  transmission  line  model  for  arc  current 
o  Parametric  representation  matching  assumed  arc  current  components  to 
measured  EM  fields  at  distant  locations 

The  physical  analog  model  describes  the  detailed  physics  of  the  lightning 
channel  in  terms  of  equations  of  conservation  of  mass,  momentum,  and  energy, 
equations  of  state,  and  Maxwell's  equations.  This  type  of  model  requires  a 
detailed  knowledge  of  physical  parameters  such  as  the  ionization  and 
recombination  coefficients  and  of  thermodynamic  properties  such  as  the  thermal 
and  electrical  conductivities.  Using  this  basic  approach,  one  can  attempt  to 
predict  the  channel  current  as  a  function  of  height  and  time.  From  a 
knowledge  of  the  current,  the  remote  electric  and  magnetic  fields  can  be 
calculated  (e.g.,  Uman  et  al..  Reference  4-1  and  4-17).  Modeling  of  this  type 
has  recently  been  attempted  for  lightning  return  strokes  by  Strawe  (Reference 
4-2)  and  Gardner  (Reference  4-3)  and  holds  considerable  promise  for  providing 
a  better  understanding  of  the  return  stroke.  At  present,  such  modeling  is 
limited  by  assumptions  necessary  to  define  the  physical  parameters. 


A  less  sophisticated  level  of  modeling  Involves  mathematically  describing  the 
lightning  channel  as  an  R-L-C  transmission  line  with  circuit  elements  that  may 
vary  with  height  and  time.  The  Intent  of  these  models  Is  to  predict  arc 
channel  current  as  a  function  of  height  and  time,  and  to  use  this  current  to 
calculate  the  EM  fields.  Price  and  Pierce  (Reference  4-4)  and  Little 
(Reference  4-5)  have  used  this  approach  for  return  strokes. 

In  the  least  sophisticated  approach  to  modeling,  and  that  which  has  been  used 
extensively  for  analysis  of  lightning  measurements,  a  temporal  and  spatial 
form  for  the  channel  current  Is  assumed  and  then  used  to  calculate  the  remote 
fields.  The  assumed  current  Is  constrained  in  Its  characteristics  by  the 
properties  of  lightning  currents  measured  at  ground  level  and  by  the  available 
data  on  the  measured  electric  and  magnetic  fields.  Lin  et  al.  (Reference 
4-6)  have  reviewed  the  literature  on  this  last  type  of  modeling  for  return 
strokes  and  have  presented  a  new  return  stroke  model  which  Is  superior  to 
previous  models  of  this  type. 

In  the  next  sections,  particular  lightning  models  of  the  three  types  are 
reviewed  and  evaluated.  A  summary  of  the  major  models  Is  given  in  Table  4.1. 

The  model  chosen  for  the  Initial  analysis  Is  Strawe's  model  because  of 
simplicity  compared  to  Gardner's  model.  It  Is  less  expensive  to  run  than 
Gardner's  and  It  Includes  physics  of  the  lightning  channel.  Also,  the  results 
to  date  Indicate  agreement  with  previous  spark  modelling  and  experimental 
measurements  of  spark  temperatures  vs.  time  and  radius  vs.  time.  More 
complete  comparisons  with  velocity  of  propagation  vs.  altitude  are  being  done 
as  well  as  comparison  with  EM  field  measurements. 

The  Strawe  model  was  chosen  for  use  over  Uman's  field  fitting  model  because  of 
the  desire  to  predict  current  profiles  at  aircraft  altitudes.  The  Uman  model 
does  not  have  this  capability. 

4.1  TRANSMISSION  LINE  MODELS 

The  many  analytical  models  for  predicting  the  arc  current  of  lightning  return 
strokes  have  Involved  linear  charged  transmission  lines.  Models  include  B  vce 
and  Golde  (Reference  4-7),  Uman  and  McLain  (Reference  4-8),  Price  and  Pierce 
(Reference  4-4),  Rosich  (Reference  4-10),  and  most  recently,  Little  (Reference 


4-5).  These  models  predict  reasonable  peak  currents,  pulse  durations,  and 
propagation  velocities  provided  the  proper  time  Independent  lightning  channel 
parameters  are  chosen.  However,  since  the  lightning  channel  expands  in 
diameter  by  over  an  order  of  magnitude  during  a  pulse,  constant  parameters 
cannot  adequately  represent  the  entire  pulse  period.  More  importantly  for  a 
threat  model,  the  linear  time  independent  line  model  produces  a  step  rise  in 
current  which  contains  no  information  on  limits  for  current  rise  rate.  Rise 
rate  determines  the  level  of  induced  transients  in  exposed  equipment  during  a 
lightning  pulse. 

The  return  stroke  of  *  lightning  discharge  in  a  cloud-to-ground  flash  can  be 
represented  in  a  very  simple  way  if  the  charge  on  the  initiating  leader 
channel  is  ignored.  The  capacitance  of  the  cloud  then  contains  all  the 
electrical  charge  present,  and  the  leader  channel  is  a  resistive  and  inductive 
element  in  a  simple  series  ICR  circuit.  The  final  step  of  the  leader  acts  as 
a  switch  which  completes  the  circuit.  The  representation  cannot,  of  course, 
provide  any  information  about  the  progress  of  the  return  stroke  current  pulse 
along  the  leader  channel.  The  current  pulse  shape  is  regarded  as  identical  at 
all  points  along  the  channel. 

The  leader  channel  may  be  treated  as  a  uniform  transmission  line  as  a  first 
step.  Price  and  Pierce  (Reference  4-4)  describe  the  development  of  this  model 
(neglecting  the  cloud  capacitance)  from  the  first  discussions  of  a  line 
excited  by  a  source  at  the  base  by  Bruce  and  Golde  (Reference  4-7).  They 
consider  finally  a  finite,  lossy,  uniform  transmission  line  terminated  in  its 
characteristic  impedence  and  initially  charged  to  a  uniform  potential.  The 
current  at  the  ground  is  found  to  rise  instantaneously  to  a  maximum  value  and 
then  decay,  exponentially  at  first.  A  similar  onset  is  predicted  for  current 
flow  at  all  points  of  the  channel,  but  the  magnitude  of  the  current  peak  falls 
as  the  height  above  ground  increases. 

The  behavior  of  the  current  at  ground  level  is  in  general  accord  with 

observation,  though  peak  values  of  dl/dt  occur  at  t  =  0,  whereas  there  is  a 

delay  before  this  peak  occurs  in  real  lightning.  At  very  long  times  the 

-1/2 

current  falls  as  t  and  this  has  some  correspondence  with  the 

intermediate  current  often  observed. 


The  fact  that  the  line  Is  lossy  is  important,  since  a  loss-free  line  produces 
a  current  in  its  terminating  resistances  that  is  constant  between  abrupt 
step-changes.  Such  a  current  waveform  is  not  typical  of  lightning  current 
pulses.  Resistance  reduces  the  importance  of  reflections,  and  any  model  of  a 
lightning  channel  must  include  resistance  if  a  realistic  current  pulse  at 
ground  level  is  to  be  reproduced.  The  effective  channel  resistance  per  meter 
is  difficult  to  assess  at  present.  It  must  be  large  enough  to  prevent 

oscillations  in  the  waveform  at  ground  level,  since  the  natural  lightning 
pulse  is  unidirectional. 

Rosich  (Reference  4-9)  examines  several  approximate  techniques  for  determining 
the  effective  transmission  line  characteristic  Impedance  of  a  lightning  return 
stroke  channel.  These  formulations  Include  (1)  the  use  of  a  coaxial  line 
model  with  center  conductor  dimension  equal  to  that  of  the  central  arc  core 
for  Inductance  calculation  and  of  the  average  positive  corona  sheath  for 
capacitance,  (2)  a  vertical  line  segment  quasi  static  (monopole)  model  by 
Berger,  (3)  a  dipole  model  with  sinusoidal  current  distribution  by 

Schelkunoff,  and  (4)  a  2  dimensional  numerical  (monopole)  solution  of 
Maxwells'  equations  by  Rosich  himself.  The  models  are  all  linear  and  with 
constant  (In  time  and  position)  parameters. 

Rosich  applies  the  models  to  a  range  of  arc  and  corona  radii  to  obtain  a  range 

of  effective  characteristic  Impedances  (ZQ).  ZQ  lies  In  a  much  more 

compact  range  than  the  radial  dimensions  since  it  Is  essentially 
logarithmically  dependent  upon  them.  The  numbers  he  obtains  are  typified  by  a 
50  to  500  ohm  resistance  in  series  with  a  few  hundred  picoFarads.  They  are 
typical  of  the  model  values  commonly  assumed  in  transmission  line  channel 
model s . 

Little  (Reference  4-5)  determines  the  transmission  line  parameters  from 
physical  arguments  considering  only  return  strokes.  To  calculate  capacitance, 
the  leader  channel  is  regarded  as  a  cylindrical  charged  conductor  of  small 
radius  extending  almost  to  the  ground  with  a  gap  between  the  end  of  the  leader 
and  the  ground.  The  capacitance  Is  determined  by  the  gap  size,  dimensions, 
shape  and  height  of  the  cloud.  The  Inductance  is  essentially  that  of  a 
coaxial  line  of  small  radius  and  Is  constant  throughout  the  pulse.  The 
resistance  is  taken  large  enough  to  make  the  current  pulse  undirectional  but 


small  enough  to  keep  some  current  oscillations.  All  parameters  are  fixed 
throughout  the  entire  lightning  pulse.  Comparison  of  Little's  computed 
current  values  are  compared  to  data  In  Table  4.2.  The  peak  current  values  are 
high  compared  to  the  data. 

Lin  et  al  (Reference  4-6)  has  tested  the  ability  of  the  Bruce-Golde  and 
transmission  line  models  to  predict  near  and  distant  electric  and  magnetic 
field  data  for  early  and  late  times.  Figure  4.1  shows  the  inability  of  either 
model  to  match  the  field  data  at  all  times.  This  Inability  to  match  the  field 
data  throughout  the  pulse  shows  the  necessity  of  having  time  varying 
parameters  to  specify  l'-’htnlng  characteristics.  Also  the  Inability  to  match 
both  near  and  far  fields  simultaneously  points  out  the  desirabilty  of 
Including  self-consistent  physics  in  defining  the  current  channel. 

4.2  CURRENT  MODELS  DERIVED  FROM  FIELD  DATA 

Modeling  attempts  using  field  data  to  determine  lightning  current 
characteristics  are  discussed  In  this  section.  The  current  Is  time  and  height 
dependent,  an  Improvement  on  the  previous  transmission  line  models  discussed 
in  the  previous  section.  Lin,  Uman  and  Standler  (Reference  4-6)  have  defined 
a  fairly  complex  current  distribution.  The  channel  current  Is  decomposed  into 
three  components  (continuing  current,  breakdown  current  and  corona  current) 
each  dominating  a  separate  time  period  and  having  a  physically  reasonable 
basis.  The  continuing  or  leader  current  is  constant  in  time  and  elevation  and 
turns  on  when  the  return  stroke  Initiation  begins.  The  breakdown  pulse 
propagates  up  the  channel  with  a  velocity  Imposed  that  Is  based  on 
experimental  values.  The  corona  current  Is  caused  by  the  radially  inward  and 
then  downward  movement  of  the  charge  stored  In  the  corona.  It  Is  Initiated  at 
each  height  as  the  breakdown  pulse  front  moves  by.  The  current  profiles  used 
are  shown  in  Figure  4.2.  The  fit  to  near  and  far  field  data  Is  quite  good  in 
that  both  close  and  distant  electric  and  magnetic  fields  are  matched.  Two 
examples  are  given  in  Figure  4.3. 

Lin  et  al  compares  the  peak  current  from  the  subsequent  stroke  to 
measurements.  Their  mean  value  of  23  kA  Is  reasonable  and  agrees  with  recent 
results  published  by  Garbagnatl  (Reference  4-10).  No  comparison  Is  available 
with  rise  time,  however,  which  Is  one  of  the  most  Important  parameters  that 
the  threat  model  will  have  to  predict. 
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Major  limitations  of  this  model  Include  (1)  non-uniqueness  of  the  current 
decomposition,  (2)  artificially  Imposed  corona  current  which  leads  to  imposing 
velocity  of  propagation  as  an  input  parameter  to  the  model,  (3)  fixed  height 
dependences  built  in  the  model  and  not  determined  by  the  physics. 

The  non-uniqueness  of  the  current  decomposition  from  the  EM  field  values  leads 
to  a  range  of  possible  current  rise  rates  and  maximum  amplitude  values.  This 
range  of  values  is  due  to  the  many  parameters  Introduced  in  the  model  and  is 
not  yet  quantified. 

The  model  assumption  of  a  constant  value  of  velocity  of  current  pulse 
propagation  up  the  channel  also  leads  to  uncertainty  In  the  prediction  of 
current  parameters.  The  velocity  values  picked  to  use  In  the  model  are 
averages  of  experimental  observations  but  the  observations  themselves  show  a 
wide  range  of  values.  Also,  the  height  variation  of  the  propagation  velocity 
is  not  included  which  leads  to  further  uncertainties. 

The  artificial  height  dependences  Included  In  the  model  must  be  verified  by 
correlation  with  experimental  data  or  the  predictions  can  not  be  used  for 
extrapolation  purposes  to  aircraft  altitudes. 

4.3  BASIC  PHYSICS  MODELS 

4.3.1  Strawe 

Strawe  (Reference  4-2)  has  developed  a  transmission  line  model  of  a 
cloud-to-ground  return  stroke  in  which  the  non-linear  breakdown  physics  is 
included.  Consistent  channel  and  line  parameters  are  determined  from  the  past 
current  time  history  using  a  shock  wave  (Braginskll  type)  model  of  the  channel 
arc.  The  model  determines  peak  currents,  decay  times,  and  scale  heights 
(lengths)  similar  to  linear  models.  It  also  provides  current  rise  times  and 
rise  rates  well  In  accord  with  measured  data.  It  Indicates  that  the  velocity 
of  propagation  of  the  current  wave  declines  with  elevation,  as  observed 
photographically,  even  when  the  line  model  Is  Initially  uniform  In  temperature 
and  channel  diameter.  Most  Importantly,  it  predicts  a  substantial  decline  in 
current  rise  rate  with  elevation  or  distance  from  the  discharge  Initiation. 


The  actual  current  and  the  channel  arc  Interact  In  a  lightning  discharge  to  a 
high  degree  especially  In  the  early  phase  when  both  the  current  and  arc  plasma 
are  building.  A  self  consistent  model  Is  needed  to  describe  the  buildup  phase 
which  determines  such  Important  current  parameters  as  maximum  rate  of  rise  (or 
rise  time)  and  propagation  velocity.  Detailed  numerical  studies  (Brode, 
Reference  4-11;  Plooster,  Reference  4-12,  4-13;  Hill,  Reference  4-14)  of  arc 
development  have  been  carried  out  In  cylindrical  symmetry  for  assumed  current 
time  histories.  These  analyses  have  established  the  time  and  spatial 
development  of  the  related  plasma  parameters,  l.e.,  temperature,  pressure, 
particle  densities,  conductivity,  arc  radius,  etc.,  for  specified  currents  but 
not  for  natural  sel*  consistent  ones. 

Close  agreement  between  the  step  excited  currents  on  vertical  conducting  wires 
over  conducting  ground  and  those  on  a  transmission  line  has  been  established 
(Reference  4-15).  This  forms  a  justification  for  the  use  of  the  economical 
and  traditional  line  model  of  the  return  stroke. 

In  Strawe's  model  It  Is  assumed  that  the  branched  structure  of  the  charged 
cloud- to-ground  leader  system  can  be  represented  for  channel  current 
calculations  as  a  network  of  transmission  line  segments.  The  current  so 
obtained  can  be  used  with  an  assumed  channel  geometry  to  calculate  the 
resultant  electromagnetic  fields. 

In  the  models  solved  to  date  only  the  equivalent  line  resistance  per  unit 
length  (R)  Is  considered  non-linear  and  time  varying.  In  principle,  since  the 
channel  radius  changes  significantly  with  time  and  the  channel  is 
geometrically  more  nearly  a  vertical  monopole  than  a  linear  transmission 
line.  The  Inductance  (L)  and  capacitance  (C)  per  unit  length  should  also  be 
time  functions.  These  Involve  logarithms  of  time  varying  terms  ana  are 
themselves  weak  time  functions.  Initially  these  time  dependences  are  not 
Included.  The  channel  tortuosity  (random  path  bending)  and  corona  sheath 
charge  storage  effects  are  Included  as  modifications  to  the  L  and  C  values  or 
equivalently  as  modification  to  the  hot  channel  (R  *  0)  velocity  factor  (Vf) 
and  characteristic  Impedance  (ZQ).  Some  of  the  charge  stored  In  the  sheath 
during  leader  formation  Is  retrievable  during  the  return  stroke.  This  is 
accounted  for  by  assuming  a  larger  effective  conductor  radius  for  capacitance 
calculation  than  the  sub-centimeter  current  carrying  core. 


The  resistance  models  developed  for  use  here  are  based  on  the  spark  channel 
model  of  Braglnskil  (Reference  4-16).  This  relatively  simple  model  assumes 
that  a  conductng  channel  has  been  established  prior  to  the  Initiation  of  the 
spark  by  prebreakdown  streamer,  and/or  leader  processes.  The  resultant  arc 
radius,  temperature,  pressure,  etc.  are  determined  from  the  spark  current  time 
history.  The  current  (1)  is  assumed  to  heat  (I^R  loss)  the  Initially 
conducting  arc  plasma  to  higher  temperatures  and  tens  of  atmospheres  of 
pressure.  This  condition  produces  a  hydrodynamic  shock  wave  In  the  air 
surrounding  the  spark  channel  resulting  In  a  rapid  channel  expansion. 
Braginskii  uses  the  strong  shock  approximation  to  simplify  the  physical 
picture  of  the  expansion  process.  This  picture  produces  an  essentially 
uniform  electrical  conductivity  (determined  from  channel  temperature  and 
pressure)  which  is  nearly  constant  In  time.  The  channel  resistance  per  unit 
length  Is  determined  from  the  conductivity  and  the  arc  radius. 

The  arc  channel  geometry  is  assumed  to  be  axi symmetric  or  locally  of 

cylindrical  symmetry.  The  channel  parameters  are  described  in  terms  of  a 

2 

deposited  energy  rate  set  by  the  local  current  time  history  (1  R  loss). 
Detailed  analyses  (References  4-11  to  4-14)  have  established  representative 
radial  contours  of  these  parameters  for  lightning-like  currents.  These 
analyses  solve,  in  cylindrical  symmetry,  the  hydrodynamic  equations  of 
continuity  of  mass,  momentum,  and  energy  (transfer)  together  with  two 
equations  or  state  relating  pressure,  temperature,  mass  density.  Internal 
energy  density,  etc.  The  model  was  developed  to  include  the  lost  and 
reabsorbed  thermal  radiation  from  the  channel,  temporal  variation  of  thermal 
and  electrical  conductivity,  and  low  pressure  momentum  transfer. 

In  comparison  with  Uman's  EM  field  current  fitting  model,  the  currents 
components  Included  in  Strawe's  model  are  Uman’s  breakdown  current  and  corona 
current.  To  date,  the  continuing  current  component  is  not  yet  Incorporated. 
This  addition  would  make  the  late  time  field  values  agree  more  with 
experimental  measurements. 

Representative  model  current  time  histories  are  shown  In  Figure  4.4.  The  most 
significant  effect  shown,  from  an  electrical  transient  point  of  view.  Is  the 
rapid  reduction  In  peak  current  rise  rates  with  elevation  or  distance  from  the 
discharge  Initiation  point.  As  expected,  the  peak  amplitude  declines  slowly 

and  the  rise  time  Increases  with  elevation. 
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Comparisons  of  Strawe's  model  with  results  given  by  Plooster  (Reference  4-12, 
Figure  3)  on  temperature  within  the  core  agree  well  both  in  magnitude  and 
trend  with  radius.  The  results  also  agree  with  experimental  results  on  core 
temperatures  given  by  Orville  (see  Reference  4-14).  Preliminary  comparisons 
of  the  height  dependence  of  the  velocity  of  propagation  up  the  channel  agree 
In  trend  with  some  of  Orville's  luminosity  measurements.  Agreement  with 
previous  results  for  both  spark  modelling  and  experimental  results  indicates 
Strawe's  model  has  Incorporated  most  of  the  important  physics  known  to  date. 

4.3.2  Gardner 

Gardner  (Reference  4-3)  uses  a  nonlinear  transmission  line  model  of  the  return 
stroke  channel.  It  Is  similar  in  structure,  physics,  and  mathematical  form  to 
that  of  Strawe  except  for  reduced  channel  complexity  (one  dimensional:  no 
branching,  bending  or  network  representations  of  ground  connections,  etc.)  and 
a  more  detailed  radial  treatment  of  arc  development.  He  uses  constant  (time 
and  position)  L  and  C  line  parameters,  although  this  could  be  generalized  to 
Include  dependence  on  local  arc  and  positive  corona  radii  as  well  as  local 
current  time  history.  The  only  time  varying  parameter  Is  the  arc  channel 
resistance  per  unit  length  R  (z,  t),  as  In  Strawe's  model. 

Like  Strawe's  model  the  only  source  of  energy  Input  to  the  arc  channel  is 
2 

joule  heating  (1  R)  due  to  the  channel  current.  The  model  represents  an 
extension  of  the  Strawe  model  In  that  the  3  region  (Braginskii  shock  model) 
radial  description  Is  replaced  by  a  radia1  grid  of  256  bins  allowing  for 
shockwave  buildup  and  decline  with  a  minimum  of  geometrical  restriction  or 
approximation.  The  3  region  Braginskii  radial  distribution  is  a  good 
approximation  for  early  arc  channel  development.  Gardner’s  approach  should  be 
capable  of  calibrating  and  Improving  upon  It  especially  for  late  time  arc 
development. 

Gardner’s  solution  Involves  the  simultaneous  solution  of  the  nonlinear 
telegrapher’s  equation  set  describing  current  propagation  on  the  conducting 
channel  and  the  nonlinear  hydrodynamic  equations  describing  the  local  pressure 
p,  temperature  T,  arc  radius  a,  arc  conductivity  a,  and  ultimately  R  (z,  t) 
In  terms  of  the  local  current  time  history.  These  are  nonlinear  partial 
differential  equations  requiring  numerical  solution.  Initial  conditions 
Involve  the  values  at  discharge  commencement  of  leader  channel  voltage,  arc 
channel  radius,  temperature,  and  pressure.  The  Initial  current  (normally 
zero)  is  also  required.  77 


Instead  of  Braglnskll's  simple  equation  of  state,  Gardner  uses  a  more  complex 
set  developed  by  Plooster  (Reference  4-13).  These  equations  Include  more 
physics  and  chemistry  than  does  the  Braglnskii  model,  although  Strawe  does 
Include  some  of  the  more  complex  physics  by  using  some  curve  fits  to 
Plooster's  model. 

Comparisons  of  Gardner's  results  to  both  Strawe  and  Plooster  agree  well  In 
magnitude  and  radial  variation  for  early  times  (arc  radii  <  few  cm).  Shock 
radii,  temperatures  and  conductivities  agree  between  Gardner  and  Strawe  for 
early  time  development  of  the  arc.  The  later  time  development  does  not  agree 
due  to  different  assumptions  made  in  the  two  models.  Gardner's  shows  the 
channel  developing  more  slowly  with  an  expanding  shock  wave  propagating  away 
from  the  channel.  Strawe's  model  continues  to  have  the  channel  radius  defined 
by  the  shock  radius  which  for  later  times  implies  a  rapidly  expanding  radius. 

These  comparisons  show  both  Gardner  and  Strawe  to  have  correct  physics  for  the 
Initial  arc  development.  Strawe's  results  are  good  for  arc  radii  out  to 
several  centimeters  while  Gardner's  results  should  be  good  beyond  that. 
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Table  4.2 

Transmission  Line  Model  Comparison  to 
Ground  Current  Data 
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t  Percentage  of  observed  values  above  computed  value:  bracketed  figures  are  estimates. 
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Figure  4.2 

Current  Distribution  for  M 


Figure  4.3 

Comparison  of  Lin,  et  al.  Model  to  Data 
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5.0  STATIC  ELECTRIFICATION 


Static  electrification  of  a  conventional  aircraft  can  occur  in  various  ways  as 
Illustrated  In  Figure  5.1.  Figure  5.1(a)  Illustrates  frictional 
electrification;  as  uncharged  precipitation  particles  strike  the  aircraft, 
they  acquire  a  positive  charge,  leaving  an  equal  and  opposite  negative  charge 
on  the  aircraft  and  raising  its  potential  to  tens  or  hundred  of  thousands  of 

volts.  Charging  occurs  both  on  the  metal  structure  of  the  aircraft  and  on 

dielectric  surfaces  such  as  the  windshield.  Dielectric  surfaces  can  thus 
become  charged  with  respect  to  the  airframe.  Engine  charging,  Illustrated  In 
Figure  5.1(b),  occurs  when  flight  vehicles  are  operated  at  low  altitudes. 
Processes  as  yc.  Incompletely  understood  occur  within  the  engine  combustion 
chamber  and  cause  a  predominantly  positive  charge  to  be  expelled  with  the 
engine  exhaust.  This  causes  an  equal  and  opposite  (negative)  charge  to  be 

Imparted  to  the  aircraft  charging  it  to  potentials  of  tens  or  hundreds  of 

thousands  of  volts.  Exogenous  charging,  illustrated  In  Figure  5.1(c),  occurs 
when  the  vehicle  flies  In  a  region  of  electric  field,  such  as  that  generated 
between  oppositely  charged  regions  of  clouds;  this  field  can  cause  discharges 
to  occur  from  the  extremities  of  the  vehicle. 

The  operational  conditions  under  which  static  electrification  can  occur  depend 
somewhat  on  the  class  of  vehicle.  Since  airplanes  encounter  severe  charging 
during  operation  in  clouds  In  horizontal  flight,  electrification  can  continue 
for  considerable  periods  of  time  on  all-weather  missions.  On  jet  aircraft 
operating  at  iow  altitude,  engine  charging  can  be  an  additional  source  of 
long-term  electrification.  Helicopters  become  charged  while  flying  through 
naturally  occurring  clouds.  In  addition,  a  hovering  helicopter  can  stir  up 
snow  or  dust  thereby  generating  its  own  cloud  of  particles  to  produce 
frictional  electrification.  Thus,  helicopters  encounter  static  problems  in 
regions  where  conventional  aircraft  do  not. 

The  charging  process  Itself  produces  virtually  no  difficulty,  but  vehicle 
voltage  and  electric  fields  can  become  so  high  after  a  period  of  time  that 
electrical  discharges  occur.  It  Is  the  discharge  of  the  accumulated  static 
electricity  that  generally  produces  the  most  harmful  effects. 


5.1  NOISE  SOURCES 

An  Important  consequence  of  static  electrification  Is  electrical  noise.  The 
various  noise  mechanisms  that  have  been  Identified  are  shown  In  Figure  5.2. 
As  the  airplane  becomes  charged,  the  electric  fields  at  the  extremities  of  the 
vehicle  become  sufficiently  high  to  cause  corona  breakdown  of  the  air.  At  the 
operating  altitude  of  airplanes,  this  breakdown  occurs  as  a  series  of  very 
short  pulses  containing  energy  In  the  radio  frequency  spectrum.  These  noise 
pulses  can  couple  into  communication,  navigation,  or  digital  circuitry  to 
produce  Interference. 

Another  source  of  noise  occurs  when  dielectric  surfaces  on  the  front  of  the 
airplane,  such  as  the  windshield  and  radome,  are  exposed  to  frictional 
charging,  as  Illustrated  In  Figure  5.2.  These  surfaces  can  be  charged  by 
impinging  particles.  Since  these  materials  are  Insulators,  the  charge  Is 
bound  at  the  place  where  It  was  deposited  and  cannot  be  discharged  until 
sufficient  electric  charge  has  accumulated  to  produce  a  streamer  (a  spark-like 
discharge)  across  the  dielectric  surface  to  the  metal  airframe.  Streamer 
discharges  are  slow  In  duration,  and  Involve  the  transport  of  charge  over  a 
large  distance.  They  therefore  produce  radio  frequency  Interference  which  can 
couple  Into  susceptible  systems  on  the  aircraft.  In  some  cases,  the 
streamerlng  on  a  square  Inch  of  surface  In  a  critical  location  Is  sufficient 
to  disable  systems. 

A  third  source  of  Interference  that  often  occurs  Inadvertently  on  airplanes  Is 
associated  with  sparking  between  unbonded  adjacent  metal  sections  of  the 
aircraft.  For  example,  consider  Figure  5.2,  which  shows  a  break  In  the  wing; 
charging  processes  on  the  airframe  will  raise  the  potential  of  the  Inboard 
section  with  respect  to  the  outboard  section  until  a  spark  occurs  In  the  gap. 
This  spark  produces  a  short  current  pulse,  which  Is  also  a  source  of  noise. 
In  flight,  the  current  required  for  corona  discharge  from  the  Isolated  wing 
tip  is  supplied  from  the  remainder  of  the  airplane. 

Finally,  slowly  varying  Induction  pulses  can  be  produced  In  antennas  by  the 
passage  of  charged  particles.  This  noise  Is  of  Importance  only  at  VLF  or  ELF 
and  does  not  pose  much  of  a  problem  to  conventional  communication  and 
navigation  equipment.  With  the  advent  of  systems  operating  at  frequencies  of 
the  order  of  10  kHz,  however.  Induction  noise  should  be  considered. 


5.2 


CHARGING  PARAMETERS 


The  Interference  problem  due  to  any  of  these  noise  sources  depends  on  the 
charging  rate  of  the  plane.  The  most  Important  process  as  measured  In  flight 
tests  Is  frictional  charging.  Examples  of  both  frictional  and  engine  charging 
are  shown  In  Figure  5.3  (Reference  5-1). 

The  precipitation  charging  current  to  a  vehicle  Is  given  by  (Reference  5-2) 

1  =  q  c  v  A  «. 

P  eff  (3) 

where 

qp  =  Charge  per  particle 
'  *  Particle  concentration 
v  ■  Aircraft  velocity 

Aeff  *  Effective  Intercepting  area  of  aircraft. 

The  various  parameters  In  the  equation  and  their  Interdependencies  have  been 
studied  analytically.  In  the  laboratory  and  In  flight,  and  are  generally 
understood  for  the  operating  regimes  of  current  aircraft.  Typical  values  of 
particle  parameters  for  an  aircraft  operating  In  the  subsonic  flight  regime 
are  given  In  Table  5.4  for  two  cloud  types. 

Table  5.4  Precipitation  Particle  Parameters  (Reference  5-3) 


Cloud  Type 

S 

pi co  Coulomb 

c 

m3 

Cirrus 

1  -  10 

2  x  104 

Thunderstorm  Anvil 

1  -  35 

5  x  104 

Laboratory  experiments  Involving  the  charging  of  projectiles  fired  through  Ice 
crystal  clouds  were  conducted  to  determine  the  relationship  between  the  charge 
acquired  and  the  Impact  velocity.  The  results  of  these  experiments  Indicate 
that  the  projectile  charge  decreases  with  Increasing  velocity  as  shown  In 
Figure  5.4  (Reference  5-4).  These  results  were  further  verified  by  flight 
tests.  It  was  noted  that  the  observed  effect  n.ight  be  caused  by  the  melting 
of  the  Ice  crystals  by  the  energy  of  the  Impact,  since  flight-test  experience 
Indicates  that  clouds  composed  of  water  droplets  tend  to  charge  an  aircraft  at 


a  much  lower  rate  than  do  clouds  containing  Ice  crystals.  Thus,  If  an  Ice 
crystal  Is  completely  melted  upon  Impact,  greatly  reduced  charging  would 
result.  Thus,  the  charging  should  follow  value  of  the  unmelted  Ice  mass  as  It 
does  In  Figure  5.4. 

The  effective  Intercepting  area,  Ae^,  as  been  found  to  be  affected  by 
aircraft  speed  as  well  as  body  shape.  The  results  of  studies  of  water  droplet 
Impingement  on  airfoils  Indicate  that  the  effective  Intercepting  area  of  a 
typical  aircraft  would  vary  with  speed  as  shown  In  Figure  5.5  (Reference 
5-2).  Combining  the  results  of  Figures  5.4  and  5.5  yields  the  curve  of  Figure 
5.6  (Reference  5-2)  which  Indicates  the  predicted  charging  current  behavior  as 
a  function  of  speed.  It  Is  noted  that  because  of  Ice-crystal  melting,  the 
charging  rate  decreases  rapidly  at  speeds  above  1500  mph.  The  maximum 
charging  current  occurs  at  about  1400  mph  and  Is  only  2.6  times  the  charging 
current  at  600  mph.  This  result  Is  highly  significant  In  that  It  Indicates 
that  precipitation  static  problems  on  highly  supersonic  aircraft  are  not 
appreciably  more  severe  than  they  are  on  subsonic  aircraft. 

5.3  CORONA 

As  the  static  charging  Increases,  the  electric  potential  of  regions  of  the 
aircraft  can  be  Increased  to  the  point  that  corona  discharge  take  place. 
These  discharges  are  in  the  form  of  a  series  of  short  pulses.  The  individual 
pulses  associated  with  these  discharges  can  be  modeled  as  (Reference  5-4) 

f(t)  -  Ae-0^  (41 

where  A  Is  the  pulse  amplitude  and  a  Is  the  pulse  decay  constant.  Both  A 
and  a  are  functions  of  atmospheric  pressure,  and  hence  of  altitude.  The 
number  of  such  pulses  per  minute,  denoted  by  v,  Is  also  a  function  of 
atmospheric  pressure.  A  good  fit  to  observed  values  of  A,  a.  and  v  can  be 
obtained  by  using: 

A  -  7.90569  x  105  p0,25 

a  *  2.777  x  IQ"2  p  (5) 

v  »  3.83767  x  103  p0,48 

where  p  Is  atmospheric  pressure,  measured  In  torrs. 


(6) 


Pressure  and  altitude  can  be  related  by 
p  .  760  «p  -  (h-l°^gA2) 

where  h  Is  altitude  given  In  kllofeet  and  pis  again  In  torrs.  The  noise 
spectrum  produced  by  v  pulses  per  second  Is  given  by  (Reference  5*4) 

r  •  *  (  *  )  1/2  iJ  ♦ 

n 


( 2  ) 


Figure  5.7  (Reference  5-5)  shows  some  characteristics  of  this  spectrum. 


Note  that  the  low  frequencies  are  a  larger  threat  for  high  altitudes.  This 
trend  with  altitude  changes  for  higher  frequencies.  Above  about  10  MHz, 
higher  noise  lewis  are  present  for  lower  altitudes. 


5.4  STREAMERS 


When  charge  Is  deposited  on  dielectric  surfaces  such  as  radomes,  windshields 
or  composite  material  structures.  It  cannot  flow  freely  to  other  parts  of  the 
aircraft  because  of  the  Insulating  character  of  these  surfaces.  If  the 
potential  between  these  surfaces  and  the  main  body  of  the  aircraft  becomes  too 
great,  a  surface  streamer  discharge  will  occur. 


The  current  flow  from  a  single  pulse  of  a  streamer  discharge  can  be 
approximated  by  (Reference  5-4) 


I(t) 

where,  for  a  typical  streamer, 

a  «  0.597 

b  «  0.403 

*  0.01  A 
max 


max 


(ae”at  +  be'et) 


a  *  1.67  x  107  Hz 
0  -  3.47  x  106  Hz 


(8) 


Clearly,  this  current  Is  many  orders  of  magnitude  smaller  than  that  due  to 
lightning. 

This  waveform  has  been  used  with  a  typical  coupling  factor  of  *  3  m”1 
for  several  streamer  lengths  to  compute  the  Induced  current  with  a  wire 
located  Immediately  below  the  streamer.  The  results  of  these  calculations  are 
shown  In  Figure  5.8  (Reference  5-3). 
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The  charge  transferred  by  a  single  streamer  discharge  Is  1  to  1.5  x  IQ-9 
Coulomb.  This  Is  roughly  the  same  as  the  charge  transfer  In  a  corona  pulse. 
The  difference  In  pulseforms  produced  by  the  two  mechanisms  results  from  the 
difference  In  the  lengths  of  the  two  discharges.  The  corona  discharge  extends 
to  only  one  tip  radius  from  the  burr  or  other  Imperfection  from  which  It 
occurs.  A  streamer  on  the  other  hand  extends  many  Inches  out  on  to  the 
dielectric.  This  long  discharge  length  causes  the  streamer  to  contain 
substantial  low  frequency  energy. 

Streamer  noise  spectra  are  given  In  Figure  5.9  (Reference  5-5)  and  a  function 
of  sample  area.  The  noise  spectra  Is  proportional  to  area  since  the  total 
charge  stored  and  thus  the  discharge  current  Is  proportional  to  area. 

5.5  COUPLING  OF  NOISE  TO  AIRCRAFT  SYSTEMS 

In  general,  noise  sources  on  an  aircraft  are  located  In  one  place,  and  the 
affected  antenna  or  system  are  located  somewhere  else.  To  calculate  the 
Interference  to  a  system  by  a  noise  source.  It  Is  necessary  to  define  the 
coupling  between  the  source  and  problem  system.  Measurements  have  been  made, 
for  example,  for  a  Boeing  707  and  a  helicopter.  Using  the  measured  coupling 
values,  noise  source  spectra  can  be  calculated  and  measured.  Examples  are 
shown  In  Figure  5.10  (Reference  5-3)  for  given  source  current  levels.  Note 
that  both  coronal  and  streamer  noise  levels  are  much  larger  than  either 
daytlome  or  nighttime  atmospheric  noise  levels.  Note  also  that  for  the 
helicopter  spectra,  the  low  frequency  streamer  noise  Is  a  factor  of  two  higher 
than  the  coronal  noise. 
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PULSES 


Figure  5.2 

Noise  Generation  Sources 


Figure  5.4 

Variation  of  Particle  Charge  with  Speed 


Figure  5.6 

Predicted  Charging  Current  for  Advanced  Aircraft 


Figure  5.7 

Measured  Corona  Noise  Spectra 
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Figure  5.8 

Typical  Current  Pulses  Induced  by  Streamer  Discharges 


Figure  5.9 

Streamer  Noise  Spectra  —  Measured  and  Calculated 
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Figure  5.10 

Measured  Coronal  and  Streamer  Noise 


6.0  PARAMETRIC  THREAT  LEVELS 


The  threat  levels  presented  to  aircraft  by  both  lightning  and  static 
electricity  are  presented  here.  The  lightning  threat  levels  are  presented  for 
both  attached  and  nearby  strokes  from  statistical  ground  based  current 
measurements.  Static  electricity  threat  levels  Include  corona  and  streamer 
time  waveforms  and  frequency  spectra.  The  threat  variation  with  altitude, 
velocity,  etc.  Is  also  presented. 

6.1  LIGHTNING 

Lightning  may  be  a  single  or  multiple  stroke  event  either  cloud  to  cloud  or 
cloud  to  ground.  Cloud  to  ground  statistics  are  used  to  parameterize  the 
ground  lightning  "hreat.  Multiple  stroke  threat  characteristics  are  shown  In 
Figure  6.1.  The  number  of  strokes  and  the  total  length  of  an  event  determines 
how  long  a  system  may  be  Inoperable.  The  magnitude  of  the  strokes  determine 
whether  a  system  can  recover. 

Current  measurements  on  single  strokes  were  accumulated  by  Berger  (Reference 
6-1)  and  Garbagnatl  (Reference  6-2).  Statistical  results  are  shown  In  Tables 
6.1  and  6.2.  The  results  show  mean  values  (50%)  and  upper  and  lower  1%  levels 
for  all  parameters  listed.  The  results  are  almost  Identical.  Both  results 
were  obtained  on  nearby  mountain  tops  so  the  correlation  Is  a  good  check  on 
the  consistency  of  the  measurements.  Debate  about  whether  this  data 
represents  “typical"  lightning  Is  still  ongoing. 

TABLE  6.1 

BERGER  (NEGATIVE  STROKES) 


99% 

50% 

1% 

MAXIMUM  RISE  RATE 
(KA/ys) 

8 

40 

200 

MAXIMUM  AMPLITUDE 
(KA) 

9 

30 

120 

RISE  TIME 
(us) 

0.11 

1.1 

9 

FALL  TIME  TO  HALF  PEAK 
fys) 

20 

80 

300 

ACTION  INTEGRAL 

2.2  x  103 

6  x  104 

1.5  x 

(A2-S) 

104 


TABLE  6.2 

GARBAGNATI  (NEGATIVE  STROKES) 

99V 

50* 

1* 

MAXIMUM  RISE  RATE 
(kA/ps) 

6 

40 

220 

MAXIMUM  AMPLITUDE 
(kA) 

9 

22 

120 

RISE  TIME 
(ps) 

0.13 

1.2 

7 

FALL  TIME  TO  HALF  PEAK 
(ps) 

20 

95 

300 

ACTION  INTEGRAL 
(A2-S) 

2.4  x  103 

6  x  104 

1.5E  x  106 

To  obtain  statistical  bounds  on  the  frequency  spectra  from  this  data,  a  double 
exponential  waveform  was  fit  to  Berger's  1%,  501,  99*  levels.  The  fit  was 
made  to  match  maximum  rise  rate,  maximum  amplitude  and  action  Integrals.  The 
parameters  chosen  and  the  resulting  values  for  the  remaining  parameters  (rise 
time  and  fall  time)  are  shown  In  Table  6.3.  The  resulting  current  waveforms 
are  plotted  In  Figure  6.2.  The  associated  spectra  are  shown  In  Figure  6.3. 
Both  the  current  levels  for  a  direct  attached  lightning  stroke  and  the 
associated  magnetic  field  values  are  given.  The  threat  of  nearby  lightning 
strokes  comes  from  the  EH  fields.  The  magnetic  field  threat  level  for  a 
nearby  stroke  (range  *  50m)  Is  shown  In  Figure  6.4  for  the  same  waveforms 
given  In  Figure  6.2. 


105 


TABLE  6.3 

DOUBLE  EXPONENTIAL  FIT  TO  BERGER’S  (1975)  STATISTICS 


PARAMETERS 

LOWER 

99% 

50% 

MEAN 

UPPER 

1% 

a  (sec*1) 

2.3  X  104 

7.9  x  103 

5.2  x  103 

6  (sec-1) 

8  x  105 

1.3  x  106 

1.6  x  106 

I0  (LA) 

10. 

30.8 

125. 

MAXIMUM  RISE  RATE 
( kA/pS ) 

8. 

40. 

200. 

MAXIMUM  AMPLITUDE 
(kA) 

9. 

30. 

122. 

ACTION  INTEGRAL 
(AZ-S) 

2.2  x  103 

6  x  104 

1.5  x  106 

RISE  TIME 
(pS) 

4.4 

3.9 

3.6 

FALL  TIME  TO  HALF  PEAK 
(pS) 

37. 

90. 

140. 

The  parametric  range  of  the  cloud  to  ground  lightning  threat  Is  established  by 
Figures  6.3  and  6.4.  These  give  the  statistical  range  of  expected  values 
between  upper  and  lower  1%  values  and  the  mean  values  for  both  attached  and 
nearby  lightning. 

6.2  STATIC  ELECTRIFICATION 

The  parametric  variation  of  corona  and  streamer  noise  threat  levels  depends  on 
a  wide  range  of  variables.  The  threat  levels  depend  on  charging  rates  which  in 
turn  depend  on  surface  area  exposed  on  an  aircraft,  speed,  altitude,  and 
weather  conditions.  The  exposure  of  an  aircraft  to  this  type  of  threat 
depends  on  the  specific  aircraft  and  Its  mission  profile.  Each  aircraft's 
threat  level  will  be  different  (References  6-3,  6-4). 

Coronal  threat  levels  are  established  by  modeling  Individual  pulses  In  a 
string  of  pulse  discharges  to  obtain  a  noise  spectrum  normalized  to  charging 
rate  (Figure  6.5)  as  a  function  of  altitude  (Reference  6-5). 

Streamer  threat  levels  depend  on  dielectric  surface  area  over  which  charge  is 
stored.  Noise  threat  levels  are  presented  as  of  function  of  surface  area  in 
Figure  6.6  (Reference  6-5). 
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Variables 


Multi-Stroke  Threat 
(Non-Damage  Effects) 


Figure  6.2 

Lightning  Current  Waveforms 
Mean  and  Upper  and  Lower  1%  Values 


Figure  6.3 

Attached  Lightning  Frequency  Spectral  Bounds 
—  Current  Amplitude  and  Magnetic  Fields  — 


Figure  6.4 

Nearby  Lightning  Magnetic  Field 
—  Frequency  Spectral  Bounds  — 


7.0  METER0L0GICAI  PHENOMENON 


Operating  an  aircraft  in  all  kind  of  weather  exposes  it  to  atmospheric 
electricity  (AE)  threats.  Probability  of  aircraft  exposure  to  AE  threats  must 
be  calculated  for  various  segments  of  a  mission  (e.g.  climb,  cruise, 
descent).  The  probability  of  exposure  coupled  with  the  probability  of  the 
severity  of  exposure  determine  the  protection  level  necessary  for  individual 
vehicles. 

Probability  of  a  lightning  strike  to  an  aircraft  is  dependent  on  many 
variables.  Included  are  such  items  as  geographic,  monthly,  daily  and  height 
variations  in  the  lightning  flash  density  per  area  per  unit  of  time.  Also 
important  is  the  aircraft's  effective  cross  section  in  attracting  lightning 
strikes.  The  probability  of  a  strike  can  be  calculated  by  using  the  type  of 
formula  as  follows 


P$  -  (F)  (N)  (Aa)  (f  diurnal ^  ^  month^  ^  altitude^  ^  geographic^  ^ 

P  =  lightning  strike  probability  per  hour  of  flight  time 

5  2 
F  =  lightning  flash  rate  density  (flashes/hr/km  ) 

N  =  Number  of  strokes  per  flash 

A  =  effective  cross  section  of  aircraft  in  attracting  lightning  (km2) 
f  month  =  monthly  dependence  of  flash  rate 

f  diurnal  =  dail*  de?endence 
f  altitude  =  altitude  variation 
f  geographic  =  geographic  effect. 

Each  factor  determining  P$  will  be  discussed  separately  in  the  sections 
below. 

7.1  LIGHTNING  FLASH  RATE  DENSITY 

The  historically  and  internationally  accepted  parameter  of  lightning  incidence 
is  the  thunderstorm  day  (i.e.  a  day  on  which  thunder  is  heard  at  a  recording 
station)  (Reference  7.1).  Isobront  maps  (e.g.  Figure  7.1)  show  areas  of 
differing  levels  of  thunderstorm  incidence.  Generally  the  highest  incidence 


of  thunderstorms  is  near  the  equator  and  lowest  near  the  pole.  However, 
Figure  7.1  of  the  United  States  shows  that  large  regional  variations  exist. 
Limitations  of  thunderstorm  day  data  are  that  neither  storm  duration  or 
severity  is  taken  into  account.  These  factors  influence  the  flash  rate  which 
is  the  parameter  of  principle  interest  in  any  threat  scenario.  Information  on 
the  frequency  of  lightning  flash  occurrence  only  exists  in  regions  where 
thunderstorm  research  programs  exist.  Relating  discharge  frequency  to 
thunderstorm  days  has  been  done  for  various  regions.  Results  are  given  in 
Table  7.1. 

Examples  of  ground  flash  rate  densities  are  given  in  Tables  7.2-4  for  various 

measurement  techniques.  The  highest  density  rate  from  any  of  these  tables  is 
o 

12  flashes/km  /yr  seen  in  Table  7.3  in  Singapore  and  in  Table  7.4  in  the 

U.S.  Recently,  Maier  and  Piotrowitz  (Reference  7.2)  report  maximum  annual 

2 

ground  flash  densities  of  18  flashes/km  in  western  Florida.  Other  regions 
of  the  U.S.  were  reported  to  have  smaller  flash  rate  densities. 

Cloud  flash  incidence  must  also  be  quantified.  Cloud  flashes  are  generally 
harder  to  see  visually  or  photographically.  Also,  a  large  number  of  flashes 
cannot  be  classified  as  cloud  or  cloud-to-ground  from  visual  observations. 
Table  7.5  shows  various  measurements  of  cloud  to  cloud-to-ground  flashes  (last 
column).  The  data  ranges  from  0.9  to  9.5.  The  scatter  is  due  largely,  to  the 
problems  cited  above.  Prentice  and  Mackerras  (Ref.  7.3)  have  developed  rough 
mean  values  for  various  latitude  ranges  taken  from  twenty-nine  sources.  Table 
7.6  shows  the  calculated  values. 


TABLE  7.6 


Latitude  Range 


Average  ratio  of  cloud  to 
cloud-to-ground  flashes 


2°-19°  5.7 

27°-37°  3.6 

43°-50°  2.9 

5  2° -69°  1.8 


The  values  have  a  wide  scatter  and  are  from  relatively  few  measurement 


stations.  More  recent  data  from  the  National  Severe  Storm  Laboratory  in 
Oklahoma  (Reference  7.9)  has  shown  ratios  as  high  as  40:1  for  cloud  vs 
cloud-to-ground.  Additional  data  on  relative  occurence  of  cloud  to 
cloud-to-ground  flashes  is  needed  to  more  accurately  define  this  ratio. 

Each  flash  is  made  up  of  a  number  of  strokes  which  occur  very  close  together 
in  time.  For  cloud-to-ground  flashes,  interstroke  time  intervals  range  from 
10  to  100  ms  (Reference  7.4).  Cloud  flashes  can  have  much  shorter  interstroke 

time  intevals  of  0.1  to  1.0  ms.  Thomson  (Reference  7.4)  has  statistically 
analyzed  ground  flashes  from  a  large  number  of  workers  worldwide.  The  maximum 
number  of  strokes  per  ground  flash  is  about  4.2  at  the  50%  level  and  up  to  16 
strokes/flash  at  the  upper  1%  level.  This  is  compared  to  an  upper  flash  rate 
for  cloud  events  seen  on  aircraft  of  10  strokes/flash  (Reference  7.5)  and 
104  strokes/flash  (Reference  7.6). 

7.2  MONTHLY  VARIATION 

Variation  of  monthly  lightning  activity  is  dependent  on  latitude.  In 
temperate  climates,  the  maximum  occurs  in  mid-summer  with  some  activity 
occasionally  in  mid-winter  (Reference  7.1).  In  subtropical  thunderstorms, 
large  activity  is  seen  only  in  summer.  Large  summer  activity  is  seen  in  Table 
7.6  of  data  taken  at  Kennedy  Space  Center  in  Florida  and  in  Table  7.7  from 
Brisbane,  Australia.  In  the  tropics  several  patterns  have  been  observed,  non 
consistent.  Since  large  variations  in  lightning  incidence  exist  year  to  year, 
long  term  averages  of  these  parameters  must  be  developed  to  obtain  reliable 
statistics  on  risk  factors. 

7.3  DIURNAL  CHANGES 

Diurnal  variations  in  lightning  flashes  are  fairly  well  established.  The 
maximum  flashing  rate  occurs  in  the  late  afternoon  and  evening  with  a  minimum 
about  10  hours  earlier.  The  peak  hours  of  activity  and  amplitude  of  variation 
vary  with  location  and  month.  Examples  taken  from  Reference  7.7  are  shown  in 
Figure  7.3.  A  ten  year  study  in  Brisbane  by  Mackerras  (Ref.  7.8)  concluded  a 
maximum  activity  ocurred  between  19Q0  and  2000  hours  with  the  hours  earlier  in 
the  winter  and  later  in  the  summer. 


7.4  ALTITUDE  STATISTICS 


Aircraft  lightning  strike  incidents  plotted  as  a  functions  of  altitude  were 
published  in  1977  (Ref  7-9).  The  strike  rate  plotted  in  Figure  7.4  for 
routine  aircraft  operations  show  a  peak  rate  near  3-4  km  altitude.  This  was 
the  only  data  available  until  recently. 

During  the  last  several  years  research  from  the  National  Severe  Storms 
Laboratory  (Reference  7.10)  has  shown  two  centers  of  lightning  activity  from 
VHF  studies  of  thunderstorms  that  are  separated  in  altitude.  The  lower 
maximum  is  about  5  km  in  altitude  which  agrees  with  the  previous  data  while 
the  upper  center  is  about  12  km.  Altitude  distributions  for  one  set  of  storms 
is  shown  in  Figure  7.5.  Major  flashes  are  defined  in  this  study  as  those  with 
greater  than  30  VHF  signals  at  each  station.  All  minor  flashes  (those  with 
less  than  30  VHF  signals)  were  thought  to  be  cloud  flashes. 

The  upper  altitude  maximum  of  lightning  activity  is  consistent  with  the 
findings  of  NASA's  F106-B  flight  program  (Reference  7.11).  Flying  at  lower 
altitudes  resulted  in  few  direct  strikes  to  the  aircraft.  Strikes  were  only 
obtained  by  flying  higher,  from  8  to  12  km  with  a  maximum  number  of  strikes 
per  thunderstorm  penetration  at  11  km  seen  in  Figure  7.6. 

The  reason  behind  NASA's  difficulty  in  getting  struck  at  lower  altitudes  while 
being  hit  often  higher  up  is  not  well  understood.  Une  strong  possibility  is 
that  at  lower  altitudes  the  aircraft  must  be  near  developing  leaders  or  run 
into  a  developed  channel.  At  upper  altitudes  the  aircraft  may  be  closer  to 
the  charge  centers  in  the  clouds  and  actually  trigger  a  direct  lightning 
strike. 


7.5  GEOGRAPHIC  EFFECTS 

Lightning  flash  rates  have  been  observed  to  vary  greatly  between  various 
regions  of  the  world.  An  example  of  the  great  variation  in  thunderstorm  days 
across  the  U.S.  is  seen  in  Figure  7.1.  High  frequencies  of  thunderstorms  are 
seen  in  Florida  and  in  Wyoming  and  Colorado.  Low  incidence  rate  is  seen  on 
the  west  coast,  for  example.  Worldwide,  the  trend  is  generally  a  decrease  in 
thunderstorm  activity  away  from  the  equator. 


Attempts  have  been  made  to  correlate  lightnings  flash  rate  parameters  with 
latitude.  Pierce  (Reference  7.12)  has  represented  the  latitudinal  variation 
of  the  proportion  of  discharges  that  go  to  ground  (p)  as 


v  °-i 

Nc 


[i  +  (x/3on 


(10) 


Where  Ng  =  number  of  discharges  to  ground,  Nc  =  number  of  cloud 
discharges.  More  recently  Prentice  and  Mackerras  (Ref  7.3)  have  developed  the 
following  empirical  relation  from  data  from  13  countries. 


=  4.11  +  2.11  COS  3  X  (0<X<60°) 


(ID 


Thomson  (Ref  7.4)  has  recently  studied  interstroke  time  intervals  and  the 
number  of  strokes  per  flash  to  discover  any  systematic  latitudinal 
variations.  His  conclusions  were  that  only  interstrokes  times  showed  any 
latitude  dependencies  but  only  at  tropical  latitudes. 


Generally  the  regional  variations  in  thunderstorm  severity  seem  to  be  much 
more  evident  than  global  trends,  at  least  for  thunderstorm  parameters  measured 
to  date. 


7.6  OTHER  PARAMETERS 


The  problem  of  aircraft  avoiding  lightning  is  dependent  on  the  knowledge  of 
all  variables  that  correlate  with  lightning  activity.  Not  all  of  these 
parameters  are  known.  Studies  by  the  NASA  F106B  and  NSSL  show  much  more  needs 
to  be  understood  to  predict  lightning  strikes  to  aircraft.  NASA's  experience 
with  the  F106B  showed  that  the  highest  strike  rate  occurred  for  temperatures 
of  -40°  to  -50°C  as  seen  In  Figure  7.6  (Ref  7.11).  They  also  occurred 
near  11  km  in  altitude  (Figure  7.7).  Strikes  tended  to  occur  for  all 
turbulence  and  precipitation  intensities  but  generally  more  strikes  were  seen 
in  low  turbulence  and  low  precipitation  areas  (see  Figure  7.8). 


NSSL  in  coordination  wun  uie  «AbA  Hub  also  compiled  parameters  reiatea  to 
direct  strikes  using  the  UHF  radar  at  Wallops  Island  (Reference  7.13).  Their 
data  also  showed  high  correlation  of  direct  strikes  to  the  F106  at  low 
temperatures  (-40°C  or  less)  (Figure  7.9),  low  turbulence  levels  (Figure 
7i0),  some  correlation  to  light  precipitation  (Figure  7J.1),  and  to  low 

lightning  flash  rates  (Figure  7.12).  This  last  correlation  is  very 
unexpected.  The  reason  for  low  numbers  of  direct  hits  with  high  flash  rates 
nearby  is  unknown.  Speculation  is  that  it  is  related  somehow  to  the 

triggering  mechanism  of  cloud-aircraft  lightning. 

Other  parameters  may  also  be  important  in  studying  direct  strikes  to 
aircraft.  One  parameter  not  mentioned  yet  is  aircraft  size,  shape  and 
material.  The  effects  of  these  parameters  is  not  known.  A  volume  effect  is 
suspected  i.e.  larger  aircraft  are  assumed  to  “attract"  lightning  from  a 

larger  volume  of  space.  Shape  factors  may  assume  a  role  in  triggered 
lightning.  Certainly  attachment  points  are  effected  by  shape.  Material 
effects  are  also  unknown.  Tests  by  Grumman  (Ref  7-14)  and  McDonnell  Douglas 
(Ref  7-15)  on  graphite  composite  have  shown  little  difference  in  attachment 
points  relative  to  aluminum.  However,  the  effect  on  lightning  strikes  is 
unknown. 

Our  knowledge  of  the  mechanism  of  lightning  occurring  within  storm  systems  is 
still  in  its  infancy  and  must  progress  much  further  before  we  can  be 

reasonably  certain  of  all  the  parameters  effecting  direct  strikes  to  aircraft. 

7.7  STRIKE  PROBABILITY  TO  AIRCRAFT 

The  probability  of  lightning  striking  an  aircraft  depends  on  the  type  of 
aircraft  and  its  mission  profile.  For  example,  a  fighter  may  climb  quickly  to 
high  altitudes,  cruise  for  several  hundred  miles,  then  descend  to  low  altitude 
near  its  destination.  The  total  strike  probability  is  governed  by  the  time 
the  aircraft  spends  at  different  altitudes  and  weather  conditions.  The 
probability  is  also  effected  by  the  aircraft  size,  the  time  of  year  at  wnich 
the  aircraft  is  flying,  the  time  of  day  and  the  regional  area  through  which  it 
is  moving.  These  parameters  were  discussed  in  sections  7.2-6. 


The  total  probability  of  a  strike  to  an  aircraft  can  be  calculated  from  the 
following  equation: 


p  Total  *  !>’  *1 
_i _ 

1*1 


(12) 


where  £  T0^ai  =  total  probability  of  direct  strike 

=  probability  of  strike  under  different  conditions 
(weather,  region,  time  of  day  etc.) 
t.  =  time  spent  under  above  conditions 

Considering  only  altitude  variations,  the  aircraft  is  considered  under  the 
cloud-to-ground  lightning  threat  below  20,000  feet  and  under  the  cloud 
lightning  threat  above  that  level.  Therefore,  equation  12  becomes  simply 


p  Total  -  Pc-9  *10.  *  Pc-c  'hi  go 


tl0„  *  'high 

where  Pc-g  =  probability  of  cloud-to-ground  lightning  strike 

Pc-c  =  probability  of  cloud-to-ground  lightning  strike 
tlow  =  time  sPent  at  1ow  a^titudj  (<20,000  ft) 
thigh  *  time  spent  at  high  altitude  (>2000  ft) 

Missiles  and  helecopters  spend  all  their  time  at  low  altitude  and  so  are  only 
effected  by  cloud-to-ground  lightning.  Fighters  and  transports  which  spend 
most  of  their  time  at  high  altitude  are  effected  mostly  by  cloud-to-cloud 
lightning.  This  is  reinterated  in  Table  7.9. 

TABLE  7.9 


Approximate  Strike  Probability 
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Table  7.2 

Frequency  of  Occurrence  of  Ground  Flashes 
(Visual/Optical  Methods) 
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With  data  from  Golde  (1961).  Ref:  from  Golde,  1977. 
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Table  7.4 

Frequency  of  Occurrence  of  Ground  Flashes 
(Power  Transmission-Line  Studies) 
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Ref:  from  Proceedings  of  4th  Annual  Workshop  on  Meterologicai  and  Environmental  Inputs  to  Aviation 
System,  March  1980. 


Table  7.8 

Annual  Variation  in  Occurrences  of  Lightning 

and  Rainfall,  Brisbane* 
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Figure  7.2 

lulative  Frequency  Distributions— 
Number  of  Strokes  per  Flash 
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Figure  7.4 

Aircraft  Lightning  Strike  Incidents 
as  a  Function  of  Altitude 
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Figure  7.6 

Thunderstorm  Penetrations  and  Lightning  Statistics* 


Figure  7.7 

Thunderstorm  Penetrations  and  Lightning  Statistics* 
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Precipita tion  Intensities* 


In  uoiint  panatratad  by 
lha  F-1068  aiiplana. 


Figure  7.9 

Temperature  Inside  Storms 
and  Probability  of  Direct  Strike 
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Figure  7.10 

Turbulence  Level  and 
Probability  of  Direct  Strike11 


Figure  7.1 1 

Rain  Intensity  and  Probability  of  Direct  Strik 
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Figure  7.12 

Lightening  Flash  Rate  and  Probability  of  Direct  Strike 


8.0  SUMMRY 


This  document  provides  definition  of  the  AEHP  program  threat  environment.  The 
threat  document  will  be  revised  periodically  to  provide  the  "best"  current 
definition  of  AEH  threats.  This  Issue  summarizes  the  initial  atmospheric 
electricity  threat  to  be  used  for  the  AEHP  program.  The  experimental  basis 
for  this  threat  was  critically  reviewed.  Comparisons  with  measured 
statistical  data  were  made  as  well  an  assessment  of  the  validity  of  the  data. 
Comments  on  on-going  lightning  data  collection  programs  are  made.  This  data 
will  Impact  the  final  threat  definition.  The  Initial  AEH  threat  definition  is 
needed  to  support  the  Phase  I  environmental  Impact  tests.  A  final  AEH  threat 
definition  will  arise  out  of  additional  studies  of  existing  data  as  well  as 
new  data  as  It  becomes  available. 

The  topics  Included  In  this  document  are  definition  and  review  of  lightning 
data  sources,  data  comparison  with  the  Initial  lightning  threat  parameters, 
lightning  model  review  and  discussion,  static  electrification  threat  levels, 
parametric  threat  levels  of  both  lightning  and  static  electricity,  and 
comparison  of  lightning  with  NEMP  and  EMI  threats.  A  discussion  on 
meterologlcal  environment  effects  will  be  added  In  the  next  Issue. 
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8.1  CONCLUSIONS 

The  Initial  definition  of  the  AEHP  lightning  threat  Is  compared  to 
ground-based  direct  current  tower  measurements.  The  validity  of  the  measured 
data  was  critically  reviewed.  To  refine  the  ground  lightning  threat,  the 
available  statistical  data  sets  (Including  Berger,  Garbagnatl,  Uman,  Eriksson, 
Clanos  and  Pierce  and  Popalansky)  must  be  critically  assessed  and 
quantitatively  as  to  accuracy  and  experimental  limitations.  The  earlier  data 
(Clanos  and  Pierce  and  Popalansky)  do  not  reflect  the  fast  rise  times  recently 
measured  and  thus  the  data  base  Is  biased  toward  longer  rise  times.  Both  the 
tower  current  data  and  the  lightning  current  values  derived  from  field 


measurements  are  subject  to  uncertainties.  These  uncertainties  need  to  be 
quantitatively  assessed  for  use  In  an  updated  current  threat  level. 

The  static  electrification  threat  Is  based  on  published  measurements  of 
streamer  and  corona  noise.  Static  electrification  threat  levels  may  change  In 
the  future  as  new  avionic/electronic  technology  (fly-by-wire)  and  materials 


m 


{graphite/epoxy,  Kevlar,  etc.)  are  used  In  newer  aircraft.  Charging  rates  and 
locations  of  corona  and  streamers  may  change  as  well  as  the  coupling  to 
aircraft  systems.  As  advanced  technology  aircraft  are  designed,  static 
electrification  threat  levels  may  have  to  be  updated. 

The  conclusions  reached  to  date  In  defining  an  appropriate  atmospheric 
electricity  threat  are: 

1.  Present  lightning  threat  Is  based  on  currently  available  Industry  draft 
standards  and  Boeing  Interpretation  of  newer  high  rate  of  rise  data. 
Static  electrification  Is  defined  from  prior  measurements  of  corona  and 
streamer  noise. 

2.  Newer  measurements  of  lightning  current  trends  toward  higher  current  rise 
rates.  Reinterpretation  of  older  tower  data,  newer  Instrumentation  and 
derivation  current  parameters  from  field  measurements  generally  point  to 
faster  current  rise  times  and  higher  current  rise  rates. 

3.  Airborne  data  Indicates  lightning  current  amplitudes  smaller  than  ground 
based  measurements  but  rates  of  rise  may  be  the  same. 

4.  Technical  problems  arise  In  assessing  old  and  new  lightning  measurements: 

a.  Experimental  limitations  difficult  to  quantify  from  published  results. 

b.  Sample  size  of  data  from  each  experimenter  Is  too  small  for  a 
reasonable  statistical  confidence  level. 

c.  Data  Is  often  not  presented  In  form  to  correlate  lightning  parameters 
of  Interest  (e.g.  amplitude  vs.  rate  of  rise). 


d.  Tower  measurements  may  underestimate  current  rise  rates  due  to 
Impedance  and  leader  effects. 


e.  Current  parameters  derived  from  field  measurements  are  subject  to  a 
variety  of  uncertainties.  Velocity  of  propagation  assumed  can 
Introduce  large  uncertainty.  A  factor  of  two  may  also  arise  from 
assuming  the  stroke  inltltes  at  the  ground  and  propagates  upward. 

8.2  RECOMMENDATIONS 

The  atmospheric  electricity  threat  directly  affects  aircraft  protection. 
Uncertainty  In  the  threat  level  imposes  a  penalty  factor  for  all  future 
technology  aircraft.  New  materials  (graphite  composite,  Kevlar,  etc), 
fly-by-wire  flight  control  and  Increased  systems  Integration  Introduce  new 
requirements  for  protection.  An  accurate  atmospheric  threat  level  imposes  the 
least  overprotect'  ■'n  requirements  with  consequent  lower  cost  and  weight 
penalties.  The  present  uncertainty  In  the  atmospheric  electricity  threat  is 
estimated  to  be  within  a  factor  of  two  from  the  present  AEH  threat 
definition.  The  only  way  to  Increase  the  accuracy  of  the  threat  Is  by 
obtaining  more  data.  This  objective  can  be  accomplished  In  the  following  ways: 

1.  Airborne  data  collection.  This  Is  the  most  critical  need  for  defining  an 
airborne  threat  level.  Flight  tests  should  be  continued  to  establish  a 
future  data  base.  This  method  Is  expensive  and  will  not  yield  a  large 
enough  volume  of  data  for  many  years.  The  Immediate  benefit  would  be 
better  understanding  of  the  Interaction  between  naturally  occurring 
lightning  and  aircraft  In  flight. 

2.  Ground  strike  current  measurements.  Continuation  of  studies  similar  to 
Berger,  Garbagnatl  and  Eriksson  on  towers.  The  towers  need  calibration  to 
establish  effects  of  the  tower  Inductance  and  local  ground  Impedance  on 
the  lightning  current  waveforms.  Geographical  effects  on  lightning  could 
be  established  by  a  network  of  towers.  Rocket  triggered  lightning 
experiments  offer  the  advantage  of  near  certain  strikes. 

3.  Simultaneous  measurements.  This  would  allow  correlation  of  fields,  visual 
pictures  of  the  stroke,  and  luminosity  measurements  as  a  function  of 
time.  These  would  help  establish  a  more  accurate  means  of  deriving 
current  parameters  from  field  data.  The  pictures  would  establish 
orientation  and  tortuosity  profiles.  The  luminosity  data  would  establish 
the  velocity  of  propagation  more  accurately. 


4.  Luminosity  measurements.  More  data  Is  needed  to  establish  velocities  of 
propagation  near  the  ground  as  a  function  of  distance  away  from  the 
initiation  point,  velocities  are  needed  to  establish  more  reliable  means 
of  obtaining  current  parameters  from  field  data  and  for  checking 
theoretical  lightning  models. 

b.  Reprocess  the  existing  data  base.  The  most  effective  way  to  establish  a 
more  accurate  AEH  threat  Is  to  compile  and  evaluate  all  existing  lightning 
data  to  date.  This  approach  could  bring  the  uncertainty  down  to  within  a 
factor  of  two  from  the  present  ALHP  threat.  To  analyze  the  presently 
available  data  more  appropriately: 

a.  Raw  data  must  be  collected  and  experimental  setups  and  limitations 
established  from  unpublished  sources. 

b.  Quantify  errors/uncertainties  to  put  the  experimental  data  on  a  common 
basis. 

c.  Derive  the  statistical  threat  levels  from  the  total  data  base 
accumulated. 

A  complete  lightning  definition  program  would  pursue  all  the  items  listed 
above  as  cost  allowed.  The  most  effective  approach,  however,  is  item  (51. 
This  approach  also  has  the  possible  advantage  of  establishing  future 
calibration  needs  for  various  types  of  experiments.  This  could  lead  to  better 
quality  data  being  obtained  from  future  efforts  to  quantify  naturally 
occurring  electricity  environments. 
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